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§1. Introduction: history of the problem, the 
one-dimensional Schrodinger operator 

Already in the 18th century (1742) Euler paid attention to the substitutions 
which transform solutions of a certain ordinary linear differential equation into 
solutions of another equation associated with it. Let 

L = -dl + u{x) (1) 

be the Sturm- Liouvillc operator defined on the real line (a;), and let A) be a 
solution of the linear equation Lij) = Xip. Let v{x) be an arbitrary solution of the 

Riccati equation 

u{x)=v^+v'^. (2) 
Under these conditions we have the following result. 

Lemma 1 The new function tp = tpx — v{x)ij} satisfies the new equation Lip = 
Xip, where 

L = —d^ + u{x), u{x) = —Vx + . 

Proof This lemma easily follows from the fact that L is in fact the product of 
two non-commuting factors of the first order: 

L = -dl^u = -{dx + v){dx - v), 

where u = + v'^ . Further, we have 

-tpxx + uip = Xip = -{dx + v){dx - v)ip = - {dx + v)ip. 

Hence, _ _ 

Lip = - {dx- v){dx + v)ip = X{dx - v)ip = Xip, 

where 

L = -{dx - v){dx +v) = -dl + u{x). 
We have proved Lemma 1. □ 



Hence, the Euler substitution is based on the representation of L as the 
product of first-order operators which are formally adjoint to each other if u, u e 
R: 

L=-{dx+v){dx-v) = QQ+, (3) 
and is given by interchanging them: 

QQ+ = L^L = Q+Q, ip^i> = Q+ip. (4) 

In the literature the Euler substitutions are often called 'Darboux trans- 
formations' or 'Backlund-Darboux transformations' because of the role they 
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play in the theory of non-hnear systems of Korteweg-de Vries (KdV) type for 
functions u{x,t) that depend on the time as the parameter. 

The classical mathematicians put it like this: the Darboux transformation 
L I— > L is determined by one chosen solution ip of the equation Lip = 0, and 

= i>x - {Vx/'P)lp- (5) 

We have 

v{x) = (fix/'P, where Vx + v"^ — u{x). (6) 

The algebraic language, that is, the correspondence of this with the factorization 
a + L = QQ^, was used after the works of physicists of the 30s and 40s (Dirac, 
Schrodinger, Infeld). 

^From the point of view of formal spectral theory, that is, of local solu- 
tions of the equation Lip = Xip, the Darboux-Euler transformation generates a 
transformation of the eigensubspaces for all A: 

B\ ip I — > -0, 

where Lip = Xip, Lip — Xip, ip = ipx — vip. 

We should talk about the Darboux-Euler transformation of the whole linear 
hull (with respect to all A) 

B ^ a^ipi I — > ^ aiipi, 

where Lipi — Xiipi and a,; are arbitrary coefficients. 

The transformation B is not, strictly speaking, an isomorphism: for A = 
and ip = ip we have Bip = Q, v = ^px/^- Hence, the kernel of this transformation 
is one-dimensional. 

Turning to the present 'global' spectral theory of the operator L in the 
Hilbert space £2(R), we see that if the function 1^9(2;) does not belong to £2(K), 
then the Darboux transformation B can generate an isomorphism of the spectral 
theories of the operators L and L; sometimes B \s & monomorphism on a part of 
the spectrum of L that covers all but one eigenfunction in £2 (R) ; sometimes (if 
^p{x) ^ and G £2(K)) B annihilates exactly one state of the spectrum t^, the 
lowest one, the 'basic' one in £2(]R)- We do not know the complete classification 
of all possible cases from the point of view of spectral theory in £2(R)- 

Definition 1 We call the Darboux-Euler transformation the latent algebraic 
symmetry of the spectral theory of the one-dimensional Schrodinger operator, 
or the discrete symmetry. 

We recall that the 'discontinuous spectral symmetries' or 'isospectral defor- 
mations' of the operator L are the well-known systems of the theory of solitons 
like KdV, which generate a gigantic commutative continuous group and in which 
KdV is a one-parameter subgroup, see 1^, [D- 

The Darboux-Euler transformation depends on a parameter, which appears 
in the factorization (3), that is, on the choice of a solution of the Riccati equa- 
tion (2). We denote these transformations by i?, 

BL^L = -{dx-v){dx+v), ^^i;={dx-v)ip. (7) 
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Moreover, there is still a one-parameter family of transformations of 'energy 
shift' 

Ta L I — > L + a, a — const. 

Let Ba = BoTa. 

In the theory of solitons, by means of the transformations Bq, we could create 
well-known non-reflecting potentials from the zero potential. By one application 
of the transformation Be we can obtain 'soliton potentials' from the free oper- 
ator Lq = —d^, starting with v = fctanhfc(a; — xq), ip — coshfc(a; — xq), a = fc^. 
By multiple application of the transformations Ba we can obtain all rapidly 
decreasing non-reflecting potentials, taking each time a non-vanishing solution 

'below the spectrum': 

UN + l{x) = Bq„ ° BaM-i O • • • O Bao (0), (8) 

Mo(a;) = 0, ui{x) — —2k'^ / cosh^ k(x — xo), ■■■■ 

In the literature on the theory of solitons (see by means of Backlund- 

Darboux transformations the so-called 'many-soliton potentials on the back- 
ground of finite-zone potentials' were obtained by applying the transformations 
Ba with suitable parameters to the known periodic and quasiperiodic finite- 
zone potentials. This class of potentials corresponds to the limits of finite-zone 
potentials for various degenerate Riemann surfaces. For a long time nobody suc- 
ceeded in establishing the relation between the class of non-degenerate periodic 
(quasiperiodic) finite-zone potentials (one-dimensional Schrodinger operators) 
and the theory of Backlund-Darboux transformations. We recall that finite-zone 
potentials correspond to orbits of the infinite-dimensional continuous commu- 
tative group of isospectral symmetries, the so-called higher analogues of KdV, 
which have finite dimension (Q, |^-|0). These operators have remarkable 
spectral and algebro-geometric properties: their coefficients are calculated by 
means of theta-functions of hyperelliptic Riemann surfaces, and they generate 
remarkable fully integrable Hamiltonian systems. 

An important idea on the connection between Backlund-Darboux transfor- 
mations and the theory of finite-zone potentials was proposed in 1986 (p^). 
In this work the author considered 'cyclic chains' of the Backlund-Darboux 
transformations of length + 1, that is, equations of the form 

u = Baj^ o o • • • o Bag (u) 

under the condition — aw-i = • ■ ■ = ao = 0. Developing very interesting 
technical arguments in studying such chains as certain integrable non-linear 
systems, in ||l2| the supposition was put forward that for N = 2k the potential 
u{x) is always a finite- zone potential (Weiss' hypothesis). This hypothesis was 
proved in [ p^ in a more general form. 

Suppose that N = 2k and the aj may be non-zero. If '^jLq oij — 0, then 
the potential u{x) in the cyclic chain of Backlund-Darboux transformations is 
finite-zone with a Riemann surface of genus no greater than N -\- 1. 
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If X]j=o OLj — OL ^ 0, then the potential u{x) has an 'oscillator-similar' 
asymptotic form 

^ 4(jv+l)2 +'^(^)' |a;|-^oo. 

If the potential u(x) is smooth and real (without singularities), then the spec- 
trum of the operator L = —d^ + u(x) is a combination of + 1 arithmetic 
progressions with the general difference a. At the same time there is a differen- 
tial operator A^v+i of order N + 1 such that [Ajv+i, i] = a^Ar+i. 

V. E. Adler made a numerical investigation and established that the equation 
of cyclicity of the chains of transformations B^j for a ^ has non-singular 
real solutions. In particular, for N = 2 this equation is transformed into the 
Painleve-IV equation. 

Operators of the cyclic chain 

Lq = Ln+1, Li\[, . . . , Li, Lq, . . . 

have factorizations of the form aj + Lj = QjQ^ . ^From the definition of 
Darboux-Euler substitutions we have: 

For the operators Lj we consider a set of 'ground states' defined by the 
equation 

Qi-iV'oj = = Ljipo^j. 

We assume that ipoj E £2(K), and all the translations aj and the coefficients 
of the operators Qj, are real. Then the sequence of 'creation operators' 
gives the whole spectrum in £2(R) for operators Lm {M > j) according to the 
formula 




Here Lm = Lm+i{n+i), n = 2k, I e Z. 

We note that the sequence of numbers 'Y^^K=j '^k is a combination of finitely 
many arithmetic progressions with common difference X]k=o '^k — ol, since the 
set of numbers Uj is periodic. The spectrum of the operator Lm has the form 

T.K=, = A„M, J = M, M - 1, M - 2, . . . . 

These are the basic results in the theory of cyclic Backlund-Darboux chains 
for the one-dimensional continuous Schrodinger operator It would not 

be superfluous to see in the literature a publication containing a full justifica- 
tion of this beautiful algebraic picture from the viewpoint of rigorous spectral 

^ For non-trivial cyclic chains of even length a solution of the problem of classification has 
not been obtained. 



5 



theory, that is, functional analysis. So far, this problem can be considered 'more 
or less solved' on the level of the requirements of reasonable (not superrigorous) 
mathematical physics. An understanding of what happens and a set of rigor- 
ous formulae have already been obtained, but the rigorous completeness of the 
picture has not been proved in general. 

Example 1 Let iV = (a cyclic chain of period 1); we have 

• • • = L_i = Lq = Li = ■ ■ ■ , 
where Lq = Q'^Q, a + Li = QQ^. Here there arises the Heisenberg algebra 

QQ+ = Q+Q + a, Q = d^ + — , (9) 

Since Q — dx + f (x), we have 

u{x) =vx+v' ^ - + i^-j . 



The eigenf unctions ipM-j,j (mentioned above) are the eigenfunctions of the 
quantum oscillator 




Let a > 0. The equation Q^q = has a solution ■00 £ £2(1^), ipo = e~"^^/*. 
All eigenfunctions V'a/ = (Q''')^^V'o, = —{dx — ax/2), belong to £2(R) and 
generate the spectrum of the operator L — —d^ + a^x'^ /4, 



Am = - + Ma. 

In this example everything is clear, but we should establish in the general 
case that the eigenfunctions belong to £2(R)- 



§2. The non-stationary one-dimensional Schrodinger 
equation 

The Darboux substitution for the non-stationary Schrodinger equation 

ii>t = ~ipxx + uip (10) 

is worthy of mention (see iQ). These substitutions have been used to construct 
some exact solutions of the KP equation. Starting from an exact solution 

iipt = -(fxx +u{x,t)ip, (11) 

we define the transformation 

B^^i; = ^Px-—^. (12) 



6 



The function ip satisfies the equation 

i-tft = -V'xx + u{x, t)^, (13) 

where u = u — 2{\.ogif)xx- 

We express this in algebraic language. 

Lemma 2 Suppose there is given a real connection of zero curvature 

Vt = dt-w, Wx = dx-v, [Vt,V,]=0 (14) 

such that (10) has the form 

iVt^ = ^V.V', = -{dx + v). (15) 

Then the Darboux transformation has the form 

^ = V.V'- (16) 

The new function i}{x,t) satisfies the equation 

iVtV^ = V,^V^. (17) 

Proof ^From (15) it follows that 

We apply to both sides of this equation and use the relation of zero curvature 
V^Vt = VtVs. We obtain 

iVo^VtV = Vx^i> = ?Vt(v^v) = iVt^. 

We have proved Lemma 2. □ 

The representation (15) is determined by one solution y> of (11), since the 
connection of zero curvature has the form 

Vt = at - ^ , V, = a, - ^ . (18) 

Hence, a substitution of Euler-Darboux type here hardly differs from the case 
of the stationary Schrodinger operator. 

We now consider cyclic chains of length 1. Let the Schrodinger equation and 
a chosen solution y> be such that in one step we come to this (gauge equivalent) 
equation: 

u{x,t) =u - 2{\0g{p)xx = u- C{t), 
iipt = -ipxx + u{x,t)ip. 

We obtain the relation 

ip{x,t) = eyipi —— + Ax + B], (19) 
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where A, B, C are arbitrary functions of t. For C ^ we obtain an equation of 
the form 

'C{x-xo{t)f 



^(a;,t)=cxp(^^ --^^\2^+D{t)y (20) 

For C = we have (A ^ 0): 

ip{x,t) = eyi^[A{t)x + B) = exp(A(t)(a; - a;o(i))). 

It is easy to prove the following lemma. 

Lemma 3 The Darboux chain for a non-stationary Schrddinger equation is 
cyclic of period 1 if and only if the potential u{x, t) has the form 

t) = {a{t)x) +/3(t).T + 7(i), 
and the chosen solution (p{x,t) is given in the form (19), where the equalities 

iC = 2{a^ -C^), 

iA = p- 2AC, (21) 
iB = -f-{A'^ + C) 

are satisfied. The dynamics of xq from (20) is as follows: 

2a^xo - P 

xo = ^ . (22) 

The proof is obtained by simple calculation. 
We consider two cases. 

Case 1. Let w + iv = C ^ 0. Then either 7^ or = 0. In the case a ^ 
we have a 'moving oscillator of variable form'. 

Lemma 4 Solutions of the form (19), where w{to) < 0, are in £2(1^) with 
respect to the variable x for any t>to if the condition 

Ima^it) < 

is satisfied. 

Proof The domain ReC < is invariant for the system (21), since for C = 
w + TO we have ti; < for tz; = by virtue of the equation 

w — 2(Ima'^ — 2vw). 

We have proved Lemma 4. □ 

Let — const > 0. Then the condition C — —a, a G R, gives a stationary 
solution for C in the domain u = ReC < 0. All trajectories in this domain are 
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periodic and the equation is easily integrated. The trajectories are given by the 

equation H(w,v) = const, where H = {{a — w)'^ + v'^)/w. 

Remark 1 A set of operators of the form Qa = dx+ax satisfies the commutation 
relations 

[Qa,Qb] = {a-b)-l, Q+ = -Q-a (23) 

for all complex a, h. The commutators of all these operators with Ha = QaQai 
where = ~Q-a, have the form 

[Qa,Ha] = {a-a)Q-a + {a + a)Qa. (24) 

The usual coherent states arc eigenvectors of the annihilation operators. The 
states we have studied are eigenvectors of the operators 

Q_ci) = -iil), ReC<0. 

Here 7 G C, and in formula (20) above 7 = Cxq. Hence, the states (20) are 
determined purely algebraically by means of the Lie algebra (23), which contains 
the oscillator Hamiltonian H = Ha- 

Remark 2 The well-known 'coherent states' of the oscillator are eigenvectors 
of the annihilation operator 

Qip^ = 'yi;^, 7 e C, 
Q = dx + ax, ^P^{x) = e-"(^-T/")'/2. 

Using these functions as the initial conditions, we obtain solutions of the 
Schrodinger equation in the particular form (20), where C = —a = const < 0. 
For Xo{t) we obtain motion along the imaginary straight line: 

xo{t) =xo{0) + 2it, 

iP^{x,t) = e-«(^-7/«-^o(t))V2+B(t)_ 

Question For t = suppose we have C(0) = q for all solutions of (20). For 
which {xo,q,q) is the set of functions ipq{x,0), where ipq{x,0) = e^^^"^"'"^ 
complete in £2(1^), where q nms along the curve from the point qo = —a to 
g = 00 in such a way that Keq < 0? How can we choose a minimal complete 
basis if this set is overfull? Naturally we can choose a curve 00 < q < —a and 
translations xo,g such that the set will be complete, q gR. 

Remark 3 We also note that for the oscillator u = o?x^ , o? = const > 0, 
we can perform the Darboux transformation, starting from the second solution 
in £2(K): 

^(a;,i)|^^o = P„(a;)e-«-'/2, a > 0, 
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where Pn{x) is an arbitrary polynomial in x of degree n (here a = C = const). 
For the potential u{x^ t) we have 

u = a'x' - 2(log = a'x' - E • 

We determine the dynamics of ip{x,t) by expanding Pn{x) in terms of the Hcr- 
mite polynomials, that is, in terms of eigenfunctions of the oscillator. There 
arises a motion of the poles of the potential u{x,t) with respect to t (or the 
motion with respect to t of the polynomial Pn{x,t), Pn{x,0) = Pn{x)) such 
that 

Pn{x,t)=J2(^~'^''ajHj{x) 
Pnix) = E ^j^j' Xj =a + 2ja. 

j<n 

Hj {x) are the Hermite polynomials defined by the formulae 

Ho = 1, Hj{x) = ((Q+)^e-«-')/e""''> 
where = dx — ax. We have 

n 

Pn{x,t) = a„e-^(2"+iM JJ(a; - xj{t)). 

i=i 

We can obtain a more extensive family of solutions of this type from initial 
conditions of the form 

(^(x,0)=P„(x)e^(°)(^-"°)'/2^ 

applying the operators Q~^q = Q-q to the initial conditions for solutions of (20). 
The dynamics of their poles can also be of interest. 

Case 2. Let C = 0. Then a = and we have 

u{x, t) = (3x + 7, ip{x, t) = ex.p(^iA{t){x — xo{t))) . 

This is a physically reasonable case when /? = E{t) and is periodic in t. In this 
case (known as the integrable case of L. V. Keldysh) we can also find a general 
solution. 

We note with interest that even in the stationary case u{x) = /3a; + 7 when 
/3, 7 are constants, it is more convenient to use the basis of non-stationary solu- 
tions (p{x, t), which have the form of plane waves for fixed time. 

Ultimately only the non-stationary Schrodinger equation is a law of nature. 
If the force Ux does not depend on time, then we can consider the Hamiltonian 
H in an appropriate gauge as stationary. As usual, in this case in place of the 
non-stationary equation 

iV't = Hip 
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we can use the Fourier method and solve the stationary problem 

Hi: = AV- 

in the space £2(K)- 

After this, we shall have the general solution in the form of finite or contin- 
uous linear combinations 

j 

However, this generally accepted approach makes sense only if the spectral prob- 
lem is sufficiently well solved: it is only a mathematical trick, nothing more. 

If the stationary problem is not sufficiently well solved, then the Fourier 
method is not appropriate. An example of such a situation is the constant 
electric field u{x) = Ex, where the basis of non-stationary solutions (plane 
waves) is much simpler, as we mentioned above, than the eigenfunctions of the 
stationary operator H (Airey functions). 

The advantage of this approach can be particularly important in the case 
when the potential has the form 

u{x) — Ex + uo{x), 

where the function uq{x) is periodic (a constant electric field is applied to the 
crystal) . Already in pq] the perspectives of the non-stationary approach in this 
case were discussed. A physically reasonable general non-stationary case is 

u{x, t) — E{x, t), 

where E{x,t) is the electric field, periodic in both variables. If the potential 
u{x,t) is doubly periodic ('zero electric current'), then we have an extensive 
family of algebro-geometric exactly soluble equations, where u{x, t) is expressed 
in terms of theta-functions of Riemann surfaces (see [0). We construct 
the Bloch solution ijj of (10): 

V(a: + ri,t) = e'Pi^i^(x,t), 

^(a:,i + To) = e'f°^°V(^,t), 

as the 'Baker-Akhiezer function' on some Riemann surface F of finite genus, 
that is, the parameters pi, po are connected by the equation of the algebraic 
curve F. These potentials were used for solutions of the 'KP equation', where t is 
renamed as y, and time, denoted by t, is the third variable. Undoubtedly these 
solutions of Krichever , ||l^ are obtained from the cyclic chains of Darboux 
transformations, but a rigorous assertion has not been proved up to now. 

More complex is the case of a non-zero 'electric flux'. The case E{t) — f] is 
well known as the integrable case of Keldysh. 

In these cases, to understand the problem it would be useful to broaden the 
number of known exact solutions. In |T5[ under certain conditions ('quantiza- 
tion of the electric flux' using an elementary cell in the {x,t) plane) Novikov 
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introduced electric Bloch functions, which arc eigenfunctions for 'electric trans- 
lations' by periods, but their analytical properties for E ^ have not been 
understood and formulated up to now. 

We transform the Schrodinger equation (10) using the gauge transformation 

^' = e^^(-'*)v, L' = e'^Le-'^, (26) 

so that 

For this we need to take Ft = u{x,t), Ai = —F^. 

Let u = f3{t)x. We have F = A{t), A = j3{t). We have a solution of form 

= e''=Vfe(t), (27) 

where 

i^k = {k- A{t))^ipk, ipk = e-'^'^^-^^^^^'dr. (28) 

For 

P{t) = E = const (29) 

we have 

ij'k = e''=^e'-^'/3(fc^-i - tf. (30) 

This is a basis in the space of solutions. 

The operators of 'electric translations', which preserve the Schrodinger equa- 
tion, 

iVt^ = (31) 

where Vt = dt+iAo, Vx = dx+iAi, are fully analogous to magnetic translations. 
For a periodic force E{x,t), 

E{x, t) = Aox - Alt = E{x + Ti,t) = E{x, t + To), 

the electric translations that commute with (31) are given by 

fi^==e^^''("'*V(2; + ri,t), fo^ = e^^i(^'*)V(a;,t + ro), (32) 

where dFo = A{t + Tq) - A{t), dFi = A{x + Ti) - A{x), A = Aodt + Ai dx. We 
have the identity 

fifo -fofi = e'*", $01=/ / E{x,t)dxdt. (33) 

Jo Jo 

If = 27r(7, g e Z, then the group of 'electric translations' is commutative (as 
in the magnetic case). In the particular case (29) we obtain 
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The electric Bloch states, which are eigenstates for Ti,To under the condition 
ETqTi = 2'Kq, q S Z, are formally written in the form 

i,{x,t,po,Pi) = e-"''P"^"fS^{^P'kix,t)). (35) 
Obviously we have 

flip = e'P^^^ip, pi = k (mod ETq). 
These are complicated generalized functions. For f3{t) — E we have, for example, 

V'(x,i,po,Pl) = e^'-^l"^ ^^i-mTo(po+E.) + 4ikE-^-t-rnTor)\ (gg^ 

(As opposed to the magnetic case, where the functions ipk rapidly decrease 
as \x\ —^ oo, and the analogue of the series (35) converges, giving analytic 
functions.) 

As the simplest example we show that under the 'strong integer' condition 
g, s, Z e Z 

EToTi = 2nq, E^/3 = 2tts, Tq = I, (37) 
we obtain for to = kE-'^, <k < ETq 

^ij{x,to,Po,Pi) = e'''''YS{To{po + Ex) + 2Trm), k = pi (mod STq). 

mGZ 

(38) 

Hence, under the condition (37) the electric Bloch states are obtained by 
solving the Cauchy problem with singular initial condition (38). In the sur- 
vey although the scheme and the ideology were correct, misprints and 
errors were made in the formulae of the text in the description of the electric 
Bloch states. 

The function (36) is a sum of two functions V' = ^+ + ^- , where -0+ is the 
boundary value of the function which is analytic for Impo > 0, and ip- is the 
boundary value of the function which is analytic for Impo < 0. This property 
is also possibly true in the general case for electric Bloch states. Our series 
(36) converges for pi — k^ + ik^ , > 0, as we can easily see from (36), so in 
the variable k — pi (mod ETq) the function (36) is the boundary value of an 
analytic function. 

In the case of a constant field under the strong integer condition (37) , starting 
from (38) the electric Bloch function can be completely calculated. We recall 
that a theta-function with zero characteristics has the form 

e^{z) = J2 e*^™'/2+'"^ (39) 
and is well defined for Imr > 0. 
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Comparing with (36) we have finally 
7p{x,t,po,pi)=9riz), T = 2{ETo)^(^^-ty + iO, 

that is, the limit of the theta function when Imr - 

Pi=k (mod ETo), 
f 17/; = e'P"^° faV' = e'P'^' 7/;, 

To e Z, e 2ttZ, ETqTi e 27rZ. 

In the remaining cases there arise generalized functions (distributions) given by 
infinite trigonometric sums that are cubic in m e Z. Apparently, quantities of 
this type have not been studied. 

Remark 4 In the case a = and f3 = const = E wc have a special family of 
solutions such that C ^ for t G R. From (21) we find that i/{2t + ip) = C, 

We have 

-a;o(0)+a;o(t) =/3^^ 

where p < 0. 

The motion of the centre xo for p = is the motion of a classical particle in 
the constant field E = (3. In this solution we have coincidence of the classical 
and quantum pictures. 

Wc have presented these elementary solutions here to demonstrate on the ex- 
ample A'^ = 1 how the problem on cyclic Darboux chains for the one-dimensional 
non-stationary Schrodinger operator should be correctly formulated. 

Definition 2 

We call the transformation 

u^u-F7(i), V^e-*s(*) (41) 

the gauge transformation of a one-dimensional non-stationary Schrodinger op- 
erator, where j{t) is an arbitrary function of time, and g = J. 

Definition 3 The chain of Darboux transformations for a one-dimensional non- 
stationary Schrodinger operator is a sequence of equations 

iipt Ljip = -'4)xx + Uj{x, t)ijj 

such that the operator Lj^i is obtained from Lj by the Darboux transformation 
and then by the gauge transformation acting on the function 7j(f): 

Uj+i{x,t) = Uj{x,t)+-fj{t), 

Ui = Ui-2{l0gipi)xa:, (42) 
iifjt = -(fjxx + Uj{x,t)(p. 



z = i{ETo)(k-Et-^), 
(40) 
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Hence, the Darboux chain is determined by a choice of solutions ipj and 
arbitrary functions of time Jjit), j G Z. In the case of a cycHc chain of period 
A'^ we have N solutions (pj and N arbitrary functions jj {t) . 



§3. One-dimensional difference operators 

We consider a difference operator L of the second order 

Llpn = C„„i?/'ti-1 + Vnlpn + C„-0n+l- (43) 

^From the theory of solitons we know that for the conservation, 

in the discrete case, of the latent algebraic symmetries of Schrodinger operators 
(such as isospectral deformations of KdV type and the Backlund-Darboux- 
Euler transformation) it is necessary to introduce the covariant translations 
expCVa;) in place of the translations T = e'xp{dx), where exp(Va;) = c„T ipn 
tpn+iCn (the 'covariant translation'). 

After this, the theory of the corresponding class of difference operators (43) 
is similar to the continuous case: there arise isospectral deformations of the type 
of the Toda chain for (43) or of the 'discrete KdV for (43) under the condition 
Vn = 0. There also arise analogues of Darboux transformations starting from 
factorization: 

a) a + L = QQ+, T+ = T-\ Q ^an + b^T, g+ = a„ + T-^bn, 

b) a + L^ Q+Q, Q = a„ + 6„r. 

For the factorization, as before, we have to solve a difference analogue of the 
Riccati equation 

Vn + + bl, a„+i6„ = c„, (44) 

Vn + a = a'^^+b'^_j^, Cn^anbn- (45) 

Factorization is also possible in the general non-self-adjoint case (direct and 
inverse) : 

L = PnT^^ +qn+ TnT = (a„ -I- bnT){Xn + UnT^^) = (c„ + dnT^^){Vn + WnT). 

The Backlund-Darboux transformations are defined as in the continuous 
case (g, H): 

B„ = BT„, B'^ = B'T^, (46) 
BL^L = Q+Q, ij^^^ Q+V, 
B' L^L' ^QQ+, yj^i^'^Qi;, 

T^Li — > L + a, (48) 
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L = QQ~^ (in all cases we first perform factorization and then permute the 
non-commutative factors, that is, the first-order operators). 

There naturally arises the problem on cyclic chains (||2l|, |2^), for example, 
on the direct transformations 



Sa„ o B^^_, o . . . o B„„ (L) = L. (49) 

It was established in |Q that if ^ aj = 0, then L is a finite-zone operator with 
Riemann surface ('spectrum') of genus g < [N/2]. If th e Y] = a ^ 0, then 
we arrive at the difference analogues of the theory [||, (TsfT 

We should take into account that the transformations and B'^ depend also 
on a continuous parameter, that is, on the choice of a solution of the Riccati 
difference equation (44) or (45), which depends on the choice of the initial 
point bo (or bo) and the set of signs (sgna„). We denote this joint parameter 
by K. The problem of an appropriate choice of signs is discussed, for example, 
in In any case, we have the relation 

B"^' B'-'L = L 



for any L, where the choice of the parameter n' is naturally consistent with k. 

In principle we can consider general words in the group generated by the 
transformations T^'B"^, i?" for all (k, a) and formulate the cyclicity condition 
with respect to them. Hence, we have the following conclusion. 

Conclusion. The set of Darboux transformations for a one-dimensional differ- 
ence operator is larger than for a continuous one because of the two different 



types of factorizations (above). In |21 the cyclicity condition was studied 
for 'direct' chains of the form (49) only. 

It is interesting to look at the case = 1, where we obtain an analogue of 
the oscillator, starting from the Hcisenberg relation: 

QQ+ = Q+Q + a. (50) 

This equation has a solution 



Q+ = l + ^/a + bnT. (51) 

The coefficients of such an operator Q cannot be real for all n. 

Lemma 5 We assume that the 'quantization condition' 

a + bn = b(no + n) (52) 

is satisfied, where no = a/6 € Z, a, 6 G M, 6 > 0. 

Then the operators Q and are well defined in the subspace of functions 
ip such that 

^P„ = 0, n + no<0, TpeC2{R). 
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This lemma is easily proved by direct verification. 

As was shown in [|2^, the ground state is defined according to Dirac: 

Q+V-o = 0. 

At the same time we have 

i:l^^l/{a^-^{k-l)\), k = n + no>Q. (53) 

Hence, V'on Poisson distribution. Applying creation operators we obtain 

the eigenfunctions 

Q^iIjq ^ Pk{n + no)V'On, 

where Pk{n) are the well-known Charlier polynomials, orthogonal with respect 
to the weight (V'on)i that is, the Poisson distribution on Z+. 

Hence, in the difference case the Heisenberg relation (50) is realized (among 
first-order operators) on the positive half- line Z+ . 

There is also another analogue, the 'q-analogue' of the oscillator, where the 
operators Q, depend on two parameters: 

Qca - 1 + ca"T. (54) 
Lemma 6 We have the relation 

Qc,aQja-l=«'(Q^i,aQc„a-l), = Ca\ (55) 

The proof consists of direct verification. 

In the following theorem was proved (parts of it were already given 
in II). 

Theorem 1 1) The transformation 

T n ^ 1 — n 

acts on the operators Qc, Q'^ in the following way: 

TQca - (56) 

2) The equation Qc'ipa — has a solution in C2{^) under the condition \a\ > 
1. The equation Q'^ipo — has a solution in Ci^) under the condition |a| < 1. 

3) The spectrum of the operator L^^a — Qc.aQt a ''^ interval < < 1 
has the form 

" ' ' (57) 

A„ = l-a2", n>\, if \a\<l. 



Analogously, the spectrum of the operator Lc^a — QtaQc.a in the interval < 
A < 1 is given by 



A„ = 1 - n > 1, if \a\ > 1, 

A„ = 1 - a^", n > 0, if \a\ < 1. 
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Problem Investigate the spectrum of the operators L = QQ'^, L = Q~^Q 
in £2(2) in the domain A > 1, if Q has the form (54). According to the 
hypothesis [E2|, this spectrum is continuous and Lebesgue (see §6 below). 



§4. The Laplace transformation and the two- 
dimensional Schrodinger operator in a magnetic 
field 

We write a general n-dimensional stationary scalar Schrodinger operator in Eu- 
clidean space in the form 

i = -^^(a„+a„(x))' + C/(a;), (59) 

where x = {x^, . . . , x"), da — d/dx". In particular, the 2-form of the 'magnetic 
field' has the form 

H =Y^ Hapdx" A dx>^, dH = 0, 

a</3 (60) 

Ha/S = daAfj — dfjAa — —Hap- 

The electric force has the form 

E^{diU,...,dnU), (61) 

that is, the potential U is defined up to a constant. By the gauge transformations 

Li — >ef Le-f, ^ ' — * e^V' (62) 

we can change the 'vector-potential': 

Aa^Aa- idaf. (63) 

In the physically comprehensible self-adjoint case we have a real magnetic 
field and the potential: Hap S R, C/ G M, and we choose a real vector-potential 
Aa G K. Then the gauge transformations have the form (62) (in the class of 
stationary operators), where / = iip(x)^ vi^) G The operator obtained by 
this transformation possesses as before the formal self-adjoint property with 
respect to the standard scalar product in the Hilbert space £2(R): 

(V'l,V'2) =^ Vi?2rf"^- (64) 

An important class consists of operators with smooth real magnetic and 
electric fields, which are periodic with respect to some lattice F C M" of rank n. 
In this case we introduce the following notions. 
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Definition 4 We say that an operator L is topologically non-trivial if tlie mag- 
netic field (60) is the non-zero cohomology class of a torus, 

7^ [iJ] e h2(T",R). (65) 



Definition 5 We say that an operator L possesses a topologically non-trivial 
electric field if the l-form of the electric field (in the non-relativistic stationary 
formalism) is the non-zero cohomology class of a torus: 

0^[E]e H\T\R), (66) 

E = Eadx", Ea = daU{x). In other words, the potential has the form 

U{x) = Eoa.x'' + Uo{x), (67) 

where the function Uo{x) is periodic with respect to the lattice, Eq = const € 
(K"\0). 

A topologically non-trivial electric field has arisen already in the one-dimensional 
case. As we showed in §2, it was appropriate to study such a situation in the 
framework of the non-stationary {x, i)-formalism, even if the electric force is 
stationary. 

For magnetic fields we use the stationary formalism. If 71 , . . . , 7„ G F is the 
basis of the lattice, then Haf}{x -I- 7j) = Ha/3{x). For the vector-potential we 
have 

Aaix -f- 7-,) = Aa{x) + da(pj{x), 

where ipj is some function of a; = {x^, . . . , a;"). 

We define 'magnetic translations' (|l^), which have been well known in the 
physical literature since the mid-60s: 

fjij{x) = + 7j)e^'^^"(^). (68) 

We have the relation 

fjffe ^e'^^'-ffef,, (69) 

where is the flux of a magnetic field through an elementary two-dimensional 
cell of the lattice in the plane (7j,7fc), or the scalar product (integral) of the 
cocycle [H] with the basis 2-cycle Zjk £ H2{T",Z). 

In the two-dimensional case n = 2 we can present the operator L in a 
factorized form, which we write as 

2L=(d + B){d + A)-h2W, (70) 

where A(z,z), B{z,z), W{z,'z) are some functions of z, z, 

z = X -\- iy, d = dx — idy. 
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The magnetic field has the form 

2H = B,- Aj. (71) 

If H, W arc real, then wc can, by using the gauge transformation, bring the 
operator to the self-adjoint form, so we have 

2L = QiQ2 + 2W, 

Qi = d + B, Q2 = d + A = -Q+, A = ~B. 

The operator L is always written in the form (70), that is, it 'admits factoriza- 
tion' in our terminology. 

The operator L also admits an opposite factorization of the form 

2L = Q2Q1 + 2V, 

where V = W + H, Q2 = d + A, Qi = d + B. 

Definition 6 The Laplace transformation is the transformation 

Li — >L, 2L = WQ2W-^Qi+2W, (72) 

where 2L = Q1Q2 + 2W. 

Lemma 7 a) If ijj is an arbitrary solution of the equation Lip = 0, then the 

function ip = Q2'>P satisfies the equation Lip = 0. 

b) If W ^ Wq = const, then for any solution tp of the equation Ltp = Xip 
the function tp = Q2'ip satisfies the equation Ltp = \ip. In this case 2L = 
Q2Q1 + 2W0. 

Proof The proof of Lemma 7 is very simple. Since Ltp = 0, we have 

-2WtP = QiQ2i> = Qi^. 

^From this it follows that 

-27P = W-^Qi^, 
-2i> = -2Q2tp = Q2W-^Qi^, 
-2Wii=WQ2W-^Qiip. 

The lemma has been proved for the case a). In the case b), for W = Wq = const, 
the proof is exactly the same. □ 

Lemma 8 The Laplace transformation for the equation Lip = commutes with 
the gauge transformation (62). If W = Wq = const, then this applies to the 
whole family of equations Lip = \tp. 

In particular, this transformation is defined on the gauge invariant quantities 
H, W: 

H = H + ^ddlogW, W = W + H. (73) 
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The proof of this lemma is obtained by direct verification. We omit it. 

Remark 5 The classicists, beginning with Laplace, considered this transforma- 
tion for the hyperbolic equation 

1) 2L = (9^ + A){dy + B) + 2W, 

2L = W{dy + B)W-\d^ +A) + 2W; 

2) 2L^{dy+B){d^ + A) + 2V, V^W + H, 

2L' = V{d^ + A)V-^{dy + B) + 2V. 

Here the pair of mutually inverse Laplace transformations acts on solutions 
of the equation Lil) = 0:_ 

1) L^L, V '-^ {dy + 

2) L^L', ^ ^ V' = (5. + 

As before, the gauge invariants have the form _ff , W or H , V, where 

2H = d^rB - dyA. 

For these invariants the Laplace transformations have the form analogous to (73): 

1) H = H + -did2logW, W = W + H; 

1 . ^ (74) 

2) H' -did2 log V, V = V + H'. 



These transformations have been studied geometrically; they have found 
applications in the theory of congruences of surfaces in (P^ )- The classicists 



of geometry (Darboux, Tzitzeica, Moutard, see ||25|, |26|) studied chains of 
Laplace transformations and carried out a series of useful formal calculations. 
They observed that the infinite chain of Laplace transformations 

• • ■ I ^-ii Lq, Li, L2, . . . , 

where Lj+i — Lj, is equivalent to a non-linear evolution system, which in the 
modern literature is called the 'two-dimensionalized Toda chain'. Intcgrability 



of it by the methods of soliton theory was discovered in |27 , and from the 
viewpoint of field theory and Lie algebras in In fact, a series of formal 

results had already been obtained at the beginning of the 20th century. 
We consider the following quantities, that is, the potentials of a chain: 

Wn = exp/„. 

^From (74) we obtain 



W„+i -Wn = e^"+^ - ef- = Hn+i, Hn+i - i/„ = My{fn). 
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By means of the substitution /„ = g„ —gn-i we obtain the 'two-dimensionaHzed 
Toda chain' 

ing„ = e^"+^-^"-e^'"-»"-\ (75) 

where n = dxdy. 

Darboux posed the problem of classifying the cyclic Laplace chains = Lq 
of period N. For iV = 2 he deduced, in particular, that solutions of the equation 

□g = sinh g (76) 

define the class of cyclic chains. 

His student Tzitzeica (@) posed the problem: let the operator Lq — Ln 
be such that Hq = Hm = 0. What kind of non-linear systems can arise in this 
case? For = 3 he showed that such cyclic chains are generated, in particular, 
by solutions of the equation 

However, in the hyperbolic case, which was considered by the classicists, there 
are nothing but formal identities. 

Well posed global problems arise only in the elliptic case dx ^ d, dy ^ d, 
which we also consider in relation to the spectral theory of the Schrodinger 
operator in electric and magnetic fields for n = 2, in particular when the fields 
are smooth non-singular real doubly periodic functions with some lattice of 
periods F on the plane R^. 

Example 2 Let = 1. We have 

Li ~ Lq, Hi = Hq + 2^/o, 

g/l ^ g/o 

Since Hi — Hq, we obtain A/o = 0. If fo is a non-singular doubly periodic 
function (that is, Wq ^ 0), then fo = const. 

Since fi = fo, we have Hi = = Ho- Thus, we obtain only a free operator 
in the zero magnetic field, but this is a consequence of the global hypothesis. 

Example 3 Let = 2. We have 

/2 = fo, H2 = Ho- 

Since Hi = Hq + |A/o, H2 = Ho = Hi + ^A/a, we obtain 

iA(/o + /i)=0, 

or 

fo = -/i + a. 
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Further, since 

g/2 ^ g/o ^ g/l 

g/i = g/o + Hi = + J?o + ^A/o, 

we obtain 

2(g/i _ e/o) = 1a/o, 



2(e-/o+«-e/'') = iA/o. 



Let g = fo — «/2. We obtain 

e^/^sinhg = -Ag. 

8 

This shows that in the non-singular doubly periodic case all the cyclic Laplace 
chains of length N = 2 are described by solutions of the equation 

76"/^ sinhfii = Ag. 

By stretching the axes {x, y) we reduce it to the sinh-Gordon equation 

Ag = sinh g. (77) 



Example 4 Let iV = 3, ilo = = Hn- Since 

Ho = H-i = H2 + ^Af2 = Ho + ^A{fo + /i + /s), 
we obtain for non-singular doubly periodic functions fj: 

/O + /l + /2 = C. 

Since 

g/o = g/3 = g/2 + = g/2 ^H^^ 1 ^(J^ + + J^) ^ g/2 + 

we come to the conclusion: ii Hq = 0, then we have 

g/3 = g/2 = g/0_ 

Hence it follows that 

2/o + /i = C, A = C2/0, 
g/2 = g/l + H^ + 1 A(/o + /i) = = + iA(/o + /i). 
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Finally, we have 

- = iA(/o + /i), 
g/o _ eC-2/o ^ 2/0) - -^A/o. 

Let fa=g + a. We obtain 

2 ^ 

If a = C — 2a, that is, a — C/3, then we finally have 
or, after stretching the scales, 

The following result was proved in 

Theorem 2 Suppose we are given a cyclic chain of Laplace transformations 

. . . , Lqi Li, L2, . . . , Lpf = 

SMc/i that the operators of the chain have smooth doubly periodic real magnetic 
field and potential. Then all the operators of the chain are topologically triv- 
ial and algebro-geometric with respect to the zero level: Ljip = 0. The Bloch 
functions of zero level with all complex quasimomenta are defined and sweep 
out a Riemann surface of finite genus] the coefficients of the operators Lj are 
expressed in terms of theta- functions of this surface {the 'Fermi- curve'), which 
is common for all operators of the chain. 

We recall that two-dimensional algebro-geometric operators were introduced 
and studied in j32) , Q ; for a review of their theory see Q . 

The idea of the proof of this theorem can be seen for = 2, where the 
requirement of cyclicity of the chain under the conditions of the theorem reduces 
to the equation Ag = sinh g. 

In studying surfaces with the topology of a torus in M.^ with constant mean 
curvature it was established that all doubly periodic non-singular solutions of 
this equation are finite-zone (algebro-geometric); see the survey |Q. Finally, 
this follows from the fact that any manifold of solutions of a non-linear elliptic 
equation on a torus (compact manifold) has finite dimension at a non-singular 
point. Further arguments, which lead to the finite-zone property of solutions of 
the equation Ag = sinhg, follow the standard scheme that was already devel- 
oped in 1974 for the KdV equation (see the survey j|]); commutative fluxes are 
linearly dependent on a finite-dimensional invariant submanifold. 
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A natural generalization of this idea for all > 2 also leads to the proof of 
the theorem, as was noted by the authors of , ||3l| . 

Definition 7 A semicyclic Laplace chain is a set of Schrodinger operators 

Lq, . . . , Li^, Ljj^i = Lj 

such that Ln = Lo + C. 

Definition 8 A quasicyclic Laplace chain is a set of Schrodinger operators 

Lo,---,L]\f, Lj^i^Lj, 

such that the extreme operators are completely factorized with a constant po- 
tential: 

2Lq — Q01Q02 + 2Vb, 2Ln — QniQn2 + 2Vn, 
Vq — const. Vn = const, 
Qj2 = d + Bj {z,z), Qji = d + Aj (z, z). 
^From the definition we have 

2L, = Qj2Qji + 2Wj, Wo = -Ho, Wn = -Hn, 

2Lj = WjQ+^W^^Qj2 + 2Wj. 

Example 5 A'' = 1 . A semicyclic chain of length 1 . We have 

2Li ~ 2Lo — 2C, Hi ^ Ho, 
ef^ = e^" + C = ef° + Hi, Hi = C. 

Hence we obtain 

A/o = 0. 

In the non-singular doubly periodic case it follows that 

fo = const, H — C. 

Hence we arrive at a Landau operator in a homogeneous constant magnetic field 
H = C and zero potential. Starting from the solutions (let C > 0) 

Qo2'ipo = 0, Lo'ipo =0, Vo = 0, 
we obtain all the Landau levels 

((3oi)?Ao = "01, 
2LiV'i = =^ 2LoV^i 2C?/'i, C = H. 

Hence, from the lowest level Qoipo = we obtain the first level ifii — QqiPo- 
These spaces of solutions are infinite-dimensional on every level, called 'Landau 
levels' for the Landau operator 

H + 2Lo = O02O01 + H, Qo2V'o = 0. (78) 
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Its levels are 

Xn = H + 2nH, i>n = {QoiT%- 



(79) 



Example 6 = 2. A scmicyclic chain of length 2. Here we have H2 = Hq. In 
the non-singular doubly periodic case, from the equality 

H2=Ho + ^A(/o + A) 

it follows that 

A(/o + /i) = =^ fo + fi=a. 

Since 

g/2 = g/o + (7 = g/l +H2 = e^' + Ho, 

g/i =g/o_^7jg_^ l^/o, 

we obtain 

g/o _ g/l =Hq-C, 

_g/i+g/o = _/fo-lA/o. 

Finally, we have 

2{ef°-ef^) = -C-^Afo, 
fi=a- fo- 

Thus we arrive at the equation 

g/o_ga-/o + ^ = _l^y^_ 

Let fo = g + a/2. We have the equality 

g+a/2 _ a/2-g . ^ = --Ao 
^2 4 

-2C- 2e"/2sinh.9 = Ag. 

This equation has an extensive family of doubly periodic non-singular real so- 
lutions. 

Solutions of the equation LjvV'at = (-^0 + C)ipN = are obtained from 
solutions of the equation LqV'o = by the chain of algebraic 'creation operators' 

tpN = Qn-1,2 ■ ■ ■ <3o,2(V'o), ioV'JV = {-C)tpN- 

Hence, the levels of spectrum e = and e = — C for the operator Lo are 
connected by these operators and are obtained from each other. It is easy to see 
that magnetic Bloch solutions are transformed into magnetic Bloch ones under 
the action of these operators if the magnetic flux is quantized (that is, the Bloch 
solutions are defined). 
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Lemma 9 If C 0, then the magnetic flux of a non-singular semicyclic chain 
with potentials that do not vanish is different from zero. 



Remark 6 Since = Hj + iA(/j) and fj is non-singular, we see that the 

fluxes of the operators Lj coincide: = [Hj\- 



Proof of the lemma Since 



1 \ 

and [e-''"] = [e^"-i] + [H], we obtain 

[g/iv] ^ [g/o] ^ [q ^ [g/o] ^KC= [e^"] + iV[i/]. 
Hence, for the flux N[H] = iV[i/o] = KG it follows that 

C^«, (80, 

where if is the area of an elementary cell. We have proved Lemma 9. □ 



(For purely cyclic non-singular chains we have C = 0, that is, the flux is 
always equal to zero.) 

The greatest interest is in quasicyclic chains, with which the main results of 
[|o), |3l| are connected. 

Theorem 3 Let the quasicyclic Laplace chain of the Schrddinger operators 

Lq, . . . , Ln, Lj^i = Lj, 
2Lq ~ Qq iQo,2, "^Ln — Qn,iQn.2 + 2CAr 

be given with doubly periodic real quantities {fj,Hj), and suppose that the mag- 
netic flux [Hj] is positive, [H] > 0. Then the operator Ln has two strongly 
degenerate discrete levels, which are isomorphic to the Landau levels and one of 
them is the lowest: 

a) LntP = CniP, QniP = 0, a = Cn, 

b) Ln^j = 0, a = 0. 

In addition, Cn = N[H]/K, where K is the area of an elementary cell. 

Here Qj2 ^ d + Bj{z,z), Qji = d + Aj{z,z). 
Proof First, solving the equation 

Qn.2iI^ = 
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according to the scheme of [Q, we obtain the eigenvalue Lati/' = CniP- If 
[H] > 0, such a solution exists and gives a huge space isomorphic to the Landau 
level according to |36 . A basis in this space, mentioned in ||], is 



the magnetic Bloch basis, which is suitable only in the case of the integer-valued 
flux [H] = 2-Km > 0; however, the result itself is true for all positive values of 
the flux. 

We obtain the second level as follows. We solve the equation Qq2'4' = 
according to 0, |H and obtain also for [H] > a huge space of functions ip 
isomorphic to the Landau level, which can be defined by using elliptic formulae 
for the integer- valued flux [H], common for and i?Ar, 

[H] = [Ho] = [Hn]. 
Next we factorize the operator Lq in the form 
2Lo ~ Qo,2Qo,i + 2i/o = 
Then we apply the creation operator N times to a solution of Qail> — 0: 

■0 = Qn-1.1 o • • • o Qo^iip. 
The new function ■0 satisfies the equation 

Li^ijj = = {Qn.,iQn.2 + Cn)'>P- 

The function ?/; belongs to £2 if belongs to £2 (K) . 

Now we find the number Cm- We have for the Laplace chain the following 
boundary conditions: 

For the mean values we obtain 

[g/.+i] = [g/.] + ^H,+^] - [e^^] + [H] = [ef«] + [j + 1)[H] = {j + 2)[H]. 
In particular, for j + 1 = iV we obtain 

ef- = [H] + [CN] = iN + l)[H], 
where [Cat] KCn or Cat = N[H]/K. □ 

Remark 7 Let the operator L be given in the self-adjoint form, 

Q = Qi=d + A, Q^ = -Q+ = d + B, 

B = -A. 

It is also convenient to choose the 'Lorentz gauge', for which all the formulae 
are written in [Q, 

Im(A^) - 0. (81) 
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Then we have 

2L = -Q+Q + 2W = -QQ+ + 2V. 



Further, if Lijj = 0, we carry out the Laplace transformation together with 
the gauge transformation 

L I — > e^Le-s, W = , 
L = -WQW-^Q+ + 2W, 

Choosing g = —f/2, we obtain the operator L = I'^Le^l'^ in the form 

I = e-f'^lef^ = -QQ+ + 2W, 

where Q = e'i' '"^Qef ''^ , Q+ = e-^/2Q+g-//2_ obtain explicitly the self- 
adjoint operator in the Lorentz gauge transformation. 

ExEimple 7 Let N = 1. Prom the definition for Lq we have 
2Lo = —QqQo = —QoQo + 2i?0) 

where 

Qo = d + A, -Q+-d + B, Ho = Wo = ef°. 

Therefore 

2L, = -(e^o/2Q+e/»/2)(e^o/2Q^g-/o/2) +2Ho = To = e-f^^'loef^/^ 
^Prom the quasicyclicity condition for A'' = 1 we have 

Hq = const, Li = Lo + 2Hq. 
We arrive at the Landau operator. 
Example 8 Let = 2. 

2Lo = -QjQo = -QoQ^ + 2Ha, Ho = e^", 

2L2 = -QtQ2 + 2C2 = -Q2Qt + 2C2 + 2H2, 
=C2+H2= e^' + H2. 

We see that e^^ = const. Further, 

C2 = e^' = +Hi= ef" +Ho + ^A/o, 
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where Hq = e-f° . Hence we have 

2C2 = -4e/« = A/o. 

Putting /o = g + a, where = C2/2, we obtain Ag = e°'{l — e^). By stretching 
the axes we bring this equation to the form 

Ag = l-e9. (82) 

The last equation has an extensive family of smooth real doubly periodic solu- 
tions. Even its solutions, dependent on one variable (for example, on x), which 
can be found by quadrature, give us a continuous family belonging to this class, 
with an arbitrary period in y: 

/^' = 2C2-4e^°, C2>0, 

U{T) = -AT + Ae\ r = /o, A = 2^2. 
The dependence of f^ix) is determined by quadrature: 

dr 

X — Xq 



x/-4e^ + At + E ' 



E> E,,,,„ ==A 1-loK- 



Equation (82) has occurred in various physical problems — in plasma theory, 
superconductivity theory, and so on. Hector de Vega informed one of the au- 
thors (S. P. Novikov) that equation (82) arose as the Bogomol'nyi reduction 
for the Ginzburg-Landau equation for the critical value of the parameter which 
separates superconductors of types I and II (see ^%). However, the application 
of this form of the Ginzburg-Landau equation at a critical point is not clear. 
In addition, the magnetic field in this model differs from ours by a non-zero 
constant, and in the theory of superconductors solutions with singularities are 
necessary. 



Hypothesis For periods iV > 5 the quasiperiodicity condition does not give 
smooth real doubly periodic solutions, apart from a constant. In any case, the 
constant (the Landau operator) does not allow any variation — it is an isolated 
solution. 

According to the hypothesis of Novikov, the resulting class of operators is 
'maximally exactly soluble' among two-dimensional Schrodinger operators with 
smooth real doubly periodic electric potential and magnetic field with a non-zero 
(let it be integer- valued) flux [H] — 27rm. 
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These operators have two discrete levels with infinite degeneration that are 
isomorphic to the Landau levels. One of those levels is basic (that is, the lowest 
one) and the second one is certainly not the nearest to the first one (at least 
one 'magnetic zone' lies between them). If there are three strongly degenerate 
discrete levels, then, as we suppose, the operator will coincide with the Landau 
operator. 

The general arrangement of the spectrum of Schrodinger operators in a mag- 
netic field with integer-valued flux, interesting aspects of the inverse spectral 
problem, and geometric and topological properties of the spectrum of these 
operators were discussed in , |Q , Q . 

If the magnetic flux is quantized, then we have the 'magnetic Bloch function' 
of the dispersion law with number m > 0, to G Z 

fiV = e'fi^iV^, f27/i = e^f-^-V, LVm =£m(Pl,P2)V'm, (83) 

= i:^{x,y,pi,p2), [ff] = $ e 27rZ, $ ^ 0. 
For definiteness, let the function ip be normalized by analogy with finite-zone 
theory §: 

The function '>pm{x,y,xo,yo,Pi,P2) under the fixed normalization {xo,yo) is 
defined on the 'phase space' 

{x,y,pi,p2) eM.^ xT\ (84) 
Moreover, the function ^ is defined up to normalization 

i^n^ f{Pl,P2)i^n- (85) 

The manifold of zeros of the function ipm, 4'm = 0, splits into two parts: 

1) 'topological zeros', where V'm(2;,y) = for fixed pi, p2', this part of the 
zeros changes under transformations (85), but the 'algebraic number of zeros' 
is the 'Chern class' c^™^ of the dispersion law with number m for the operator 
L of general position, where the dispersion laws do not intersect, 

emi(Pl,P2) em2(Pl,P2), mi 1712; 

2) 'geometric zeros', where locally at a general point we can assume that 

{Pl,P2) = 7ni{x,y), TO>0. 

As TO ^ c» the dispersion laws converge to constants: 

em(pi,P2) const = (2to + 1)[H], 

analogous to the forbidden zones of the one-dimensional periodic Schrodinger 
operator. 

Problem Find the analogue of the 'Dubrovin equations' for the set {jm{x, y)} 
for all to; solve the inverse problem of reconstructing the operator L, starting 
from the set of geometric data, in the phase space (84), analogous to the one- 
dimensional data, without entering into the complex domain of quasimomenta. 
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§5. Difference equations and Laplace transforma- 
tions. The hyperbolic case 

The simplest hyperbolic difference equation on the standard quadratic lattice F 
in R2 is 

anipn + bn1pn+Ti + Cn^n+Ts + d„^/'„+Ti+T2 = 0, (86) 

where n = (77.1,77.2) is the number of a lattice point, rij e Z, Ti = (IjO), 
T2 — (0, 1) are the basis vectors of the lattice. We also denote by Ti, T2 the 
basis operators of translations which act on functions defined on Z^: 

Tiip{ni,n2) = -0(771 + 1, f^2), 
T2ip{ni,n2) = ■4!{ni,n2 + 1). 

We can write (86) in the following way: 

Lij = (a„ + bnTi + CnT2 + dnTiT2)ij = 0. (87) 

It was already observed in |4C| ] that for an equation of this form there is 
a class of algebro-geometric exactly soluble cases with effective theta-function 
formulae. 

The operator L is defined up to multiplication by a function: L > fnL. 
We also have gauge transformations, and so together we have a system of the 
following equivalences of (87): 

L fnLgn, ■0n '-^ 9n'^'>Pn, (88) 

where /„, g„ are non- vanishing functions on the lattice. 

Lemma 10 The operator L admits a unique representation in the form 

a) L ^ fn{{l + U„Ti){l + VnT2) + Wn) (89) 

or in the form 

b) L = /;((i + 7;;r2)(i + <ri) + 0. (90) 

The coefficients Un, Vn, Wn, fn for the factorization a) are determined from the 
formulae 

a-n = /ri(l + Wn), /„M„ = bn, fnVn = C„, fnUnVn+Ti = d„ . (91) 

The proof is based on direct verification. 

Using both factorizations a) and b), we define two Laplace transformations 
of (87) LtP = 0. 

Definition 9 A Laplace transformation of the first type is a transformation 

L I > L = f„{wn{l + 7;„T2W,7^)(1 + 7i„ri) + Wn), 

^ I > tp ^ {1+ VnT2)ll). 
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A Laplace transformation of the second type is a transformation 



Lemma 11 1 The Laplace transformations of the first and second types are 
mutually inverse. 

The proof is obtained by simple substitution, as in the continuous case. 

We denote a Laplace transformation of the first type by A'^2~ of the 
second type by . Lemma 11 claims that Aj^2^Aj|+ = 1. 

It is appropriate to note that in a completely analogous way we can define 
the Laplace transformation of (86), which corresponds to an arbitrary pair of 
orthonormal basis translations 

(Ti,T2), (r2,Ti) A++,A++, 
(Tf\T2), {T2,T,-') Ar2+,A+-, 
(Ti,T2-i), {Ti\T,) A+-,A2-+, 
iT,-\T,-% {T,-\T,-') Ar2-,A2T. 
Thus, we have eight Laplace transformations Afj , Al^ , e, a = ±, and 

^'-12^*-21 — -L- 

Hence, the Laplace transformations here form a group with four generating 
elements AI2, the structure of which we do not know so far. 

Now we return to the transformation A^2^, described above. The Laplace 
transformations are well defined on the equivalence classes (88). Therefore, as 
in the continuous case, they should be registered on the invariants of the equiv- 
alence classes, which are analogous to the magnetic field and electric potential. 

Lemma 12 2 The complete set of invariants with respect to the equivalence 
(88) is given by the following basis functions: 

bnCn+Ti J. _ Cnbn+T2 f^.. 
i^ln — -. , -fi2n — 3 • (94) 

anOn+Ti Unan+T2 

All other invariants are expressed in terms of Kin, and their translations 
on the lattice. 

Proof The proof of the invariance of these quantities is straightforward: it is 

a direct elementary verification. It is convenient to carry out the proof of their 
completeness for the factorized form (89). Expressing a„,6„,c„,d„ in terms of 

fn,Un,Vn,Wn, WC obtain 

Kin = , K2n = TTi;;;^(^"""+^^""+^."")- 
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It is clear that the potential Wn is gauge invariant, and the function /„ can be 
shortened in equivalence classes. Hence, everything reduces to gauge equivalence 
for /„ = 1. Here the 'potential' Wn and the 'curvature' 

H = V„U„ Un+T^Vn+Ti = (95) 

are gauge invariant. We need to find all gauge invariants of the operator 

Lo = {1 + UnTi){l + VnT2) 

with respect to the transformations 

Lo I > T„(l + UnTi){l + VnT2)T^^ = ( ^ + Ti ) ( 1 + 

\ Tn+Ti / V Tn+T2 

This problem is easily reduced to the classification of gauge invariant mono- 
mials, that is, products of the quantities Un+kTi+iT2j Vn+kTi+iT2 for different 
(finitely many) k, I. Elementary combinatorics leads to the conclusion that all 
these monomials are products of monomials (94) and their translations in some 
powers. We omit these details. □ 

We use the basis of invariants Wn, i?n, by analogy with the continuous case. 
Lemma 13 3 The Laplace transformation for invariants has the form 

1 + Wn+Ti = (1 + Wn+T2) , 

Wn+TiWn+T2 , 

~ 1+W„+T2 ^ ^ 

tin — ~ ■ 

1 + Wn+T2 

The proof is obtained by substituting the expressions for Ki, K2 in terms of 
a, b, c, d. 

We now consider an infinite chain of Laplace transformations for the quan- 
tities 

fr(fe) _ rr(fc+l) ~(k) _ ,,,(fc+l) 

where n — (71,1,712), fc G Z, and hIi^^\ w^'n^^^ are obtained from the previous 

(k) 

ones by a Laplace transformation. Eliminating from the equations, we 

obtain a 'completely discrete analogue of the two-dimensionalized Toda chain': 

1 I (fc+s) 1 , (fc+i) (fe) (fc) 



1 , (fe+1) 1 , (fe) (fe) (fc) 

Problem Find a comparison between this completely discrete analogue of the 
two-dimensionalized Toda chain and the systems studied in |^, starting from 
the theory of the 'Yang-Baxter' equation. 
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Example 9 We consider a cyclic chain of length N = 2. After simple calcula- 
tions we obtain 

Wn+Ti+T2 = Wn^{C + W„+Ti)(C + W„+T2)(l + W„+Ti)~^(l + Wn+Ta)"^ 

Finally, we note that there is, as we showed above, the group of all hyperbolic 
Laplace transformations, generated by four basis transformations Af^, e, cr = ±, 
which corresponds to the choice of the bases (Tf , T^). In principle, to each word 
in this group corresponds the condition for cyclicity of a chain. However, the 
structure of this group is not known so far. 

§6. The Laplace transformation for elliptic two- 
dimensional difference operators on a regular lat- 
tice. Equations of a triangle, curvature 

It appears that in the elliptic case Laplace transformations cannot be defined 
for difference operators of the second order on a square lattice. To give their 
correct definition we should follow two principles. 

I. In order that continuous and discrete 'spectral symmetries' of continuous 
Schrodinger operators be conserved in the difference case, all translations should 
be considered as covariant, as in the one-dimensional case (see §3). 

II. In the two-dimensional case for elliptic operators we should replace the 
square lattice by an equilateral triangular lattice, where the following transla- 
tions have the same length: 

T±i, {TiT2)^\ (97) 

We denote an arbitrary pair of neighbouring translations (97) (which go anti- 
clockwise) by T*, Tj* and ascribe to this pair a Laplace transformation generated 
by factorization 

L = {Xn + UnTi + Z„r2*)(x„ -|- l/n-T^T^^^ + Zn-^T'T^^^) + W„, (98) 

where the self-adjoint operator L has the form 

L = an+bnTi+CnT2+dn^T2TlT2^ +bn-TxT^^ +Cn-T2T2^ +dn^TtT2T^^ , (99) 

all the coefficients are real, and bn,Cn, dn G K+- It is easy to see that the coeffi- 
cient a„ is a numerical function on Z^, while the coefficients 6„, c„, c?„ are a con- 
nection sitting on the corresponding edges of the lattice. As a result the operator 
contains an interaction of each vertex with all six nearest neighbours, and all 
edges have the same length. We start from the following observation (|^^, ji^). 
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Lemma \A ^ A real self-adjoint operator of the form (99) with non-zero coef- 
ficients bn, On, dn ulways admits a unique representation in the form 

L = QQ^ + Wn, Q = Xn + VnTl + ZnT2, (100) 

where = . T2 = ■ {AB)^ = B^A^. Moreover we have the equalities 

XnVn-Ti = bn-Ti , 

XnZn-T2 = C„_T2, (101) 
yn-TiZn-T2 = C^ri-Ti-Tj- 

In all, following , |Q , we obtain six different factorizations, which corre- 
spond to the following six pairs of periods {T^,T^): 1) (Ti,T2), 2) (r2,Tf ^Ts), 
3) iT,-'T2,T,-'), 4) (Tf \T2-i), 5) (T^^T^'Ti), 6) {T^^T^.T^), 

L = (5jQ+ + Wjn, 

j ~ 1,2,3, 4, 5, 6. j is determined modulo 6. 

Definition 10 A Laplace transformation Pj of type j is as follows: 

J J (^Q2) 

It is defined if Wjn > 0. 



Lemma 15 5 The Laplace transformations Pj and Pj+s are mutually inverse. 
The others are connected by the following relations: 

Win = W3n = W5n, 

Qi — QzTi ~ Q5T2, 



V Wn-Ti V Wn-Ti V V ^n-T2 

(103) 

In the equivalence class 

L ^ fnLfn (104) 

i/ie operators L^^\ L^^\ L^^^ are unitarily adjoint. 
For j 2, 4, 6 analogous formulae hold: 



Remark 8 In the definition of Laplace transformations we can discard the 
requirement w„ > 0, replacing it by the condition Wn 7^ 0, and set 

ZO) = Q+w-^Q, + 1, Q+ip. 
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We obtain the same transformations up to equivalence (104), but Pj and Pj+s 
are not exactly mutually inverse. We note that in the discrete case, as opposed 
to the continuous one, the condition w ^ does not mean that w has constant 
sign. 

Remark 9 The relations (103) were not noted in ji^-Q. 

Remark 10 If the potential Wn = wq is constant (for example, for the first 
factorization), then the Laplace transformation 

PiL^ Z^i) 

acts on the eigenfunctions of all levels Lip = Xip: 

Pii/j = Q+^, PiL ^ Q+Qi + wo = Z(i). 

In the framework of the equivalence class (104) of real self-adjoint operators 
we can always arrange that Wn = const if Wn = exp /„ ^ 0. However, these 
transformations realize a formal equivalence of the spectral theories only on the 
zero level Lip = 0, and therefore we can apply the Laplace transformation to all 
levels, generally speaking, only once (for special cases see below). 

The proof of the lemma can be obtained by direct verification. 
An investigation of the cyclic chains of difference operators has not been 
carried out. The classes of purely factorized operators of the form a) and b): 

a) L = Q+Q; b) L ^ QQ+ 

are of interest, where 

Q = l + ce'i(")ri + de'^(")T2, (105) 

and 

lj{n) = Ijiui + Ij2n2 

are linear forms in n = (^1,712). 

Following iQ, we consider the 'equation of black triangles', which have 
the form {n, n + Ti , n + T2) : 

Qi^ = 0, 

and the 'equation of white triangles', which have the form (n, n — Ti,n ~ T2): 

Q+iP = 0. 

Theorem 4 1) The equation of black triangles Qip — certainly has an infinite- 
dimensional space of solutions ip G >C2(^^) if one of the following conditions is 
satisfied: 

a) lii,l22 > 0, I11I22 - II2 > 0; 
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b) ln,l22 >0, /11/22-/21 >0; 

c') III > 0, ^11/22 — ^12 > ^11(^21 — ^12); 
C") ^22 > 0, I11I22 - III > h2{l21 - /l2). 

2) The equation of white triangles Q'^'ip = certainly has an infinite- dimensional 

space of solutions tp G £2(2^) if one of the following conditions, which are ob- 
tained from the previous ones by the transformation lij —Uj, is satisfied: 

a) Zii,Z22 < 0, 111/22 - II2 > 0; 

b) /ii,Z22 < 0, ;iiZ22 - I21 > 0; 

c') ^11 < 0, ^11^22 - II2 > hi{li2 - I21); 

c") ^22 < 0, ZiiZ22 - I21 > h2{h2 - hi)- 

Proof Wc construct explicitly the solutions of our equations if the conditions 
presented above are satisfied. We examine first the case of black triangles. We 
are looking for a solution of the form 

Yn — ^ An 5 

where 

K2{n) = an\ + 2j3nin2 + 5n\. 
For Xn we obtain the equation 

(106) 

We consider three cases: 

a) the coefficients (106) depend only on m; 

b) the coefficients (106) depend only on n2; 

c) the coefficients (106) depend only on m + n2- 
This leads, respectively, to the conditions: 

a) li2 = 2p, I22 = 25; 

b) Zii = 2a, hi = 2p; 

c) hi -2a = li2 - 2/3, Z21 - 2/3 = ^22 - 26. 
In the case a) we make the substitution 

X„ = i«">„i (107) 

and obtain the equation for ipm- 

= ^m + ce-«e('"-2«)('»+i)^„+i + wde-^e^^'^-^''^"'ipm. 
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In addition we put In = 2a. After this we obtain 

(1 + «;ce-'=^/2e('=i-'i^)'")v?„ = -ce-"^m+i- (108) 
Two cases are possible. 

Case 1. hi > h2- We choose the value of Wq such that 

1 + Wgde-^/^e^^^'-^^^i = 0, (109) 

where 26 = I22, 2(3 = Z12. ^Prom this we obtain the 'quantization condition', 
which selects a discrete series of admissible values w = Wg, where q runs over Z. 

Under this condition ipm^ = ioi m > q. 
Further, we consider the relation 

= -(ce-")(l + w;,e'==-V--2'')'")-Vm+i 

as m ^ 00. If I21 — 2/3 > 0, then for m < and \m\ ^ 00 we have the 
asymptotics 

(^(;?) ~ (const) l-l. 
Our solution is constructed in the form 

=e-^2(")u;"V(9), (110) 

where 2K2{n) = hinl + 2li2n\n2 + 122^% This form is strictly positive if and 
only if 

Zii > 0, I22 > 0, I11I22 - 1I2 > 0- 

The function "^n^ certainly belongs to the space £2(Z^) and takes the value zero 
for 111 > q. 

Case 2. Z12 > hi- In this case a solution (pm of equation (108) grows no faster 
than (const) I '"I. rn ±00, for any w. Thus wc have constructed the necessary 
solution for the case' a) of black triangles. 

In the case b) we interchange the position of ui and n2- In the case c) we 
require that (f depends only on ni +n2- For example, in case c') we put 

2K2{n) = hml + 2li2nin2 + (^22 - hi + h2)n\ 

and make the substitution 

Xn = w"Vni+n2- 

For the function ipn we obtain the following equation: 

In the cases a)-c) for white triangles we proceed completely analogously. □ 

Other cases, when the solutions lie in £2(^^)5 have not been found so far. 
The form K2{n) can be written as follows. 
The case of black triangles: 
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a) 2K2 

b) 2K2 
c') 2K2 
c") 2K2 



11 


^12^ 


12 


^22^ 


11 


^21^ 


21 


^22^ 


11 




12 


^22 


11 


-/2I 



/l2 
/2I 



tl2 



21 



^12 ^21 
^22 



The case of white triangles: 

a) 2K2 

b) 2K2 
c') 2K2 
c") 2K2 



'-hi 


-h2 


-h2 


— I22 


'-hi 


-hi 


-I2I 


—h2 


'-hi 




-h2 


—h2 



-h2 

hi 

—hi + hi — h2 
-hi 



h2 
-hi 

-h2 



Remark 11 We can impose analogous conditions on the coefficients hj and 
their relation with a, (3, (5, demanding that the coefficients of equation (106) for 
Xn be dependent only on a combination of the form (ktii + rn2). If (ktii + Tn2) 
and {uni + vn2) are the basis of the lattice, then the condition of dependency 
on only one variable has the form: 



{hi - 2a)T = {h2 - 2/?)k, 
{hi - 2/3)t = {h2 - 26)k. 

Now we can seek a solution in the form 



An — ^nni+rn2' 



(111) 



(112) 



Further, for the functions (ps, s G Z, we obtain a difference equation of the form 

= ip,+ A{s)ip,+^ + B{s)ip,+r. (113) 

Only in the three cases mentioned above can this equation be easily solved: 

a) K = 1, T = 0; b) K = 0, r = 1; c) k = t = 1. 

In the other cases we obtain difference equations of order 2 and higher. We do 
not know how we could find solutions of them that lead to functions ip G £2(^^)- 



^From the theorem we have the following corollary. 
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Corollary 1 Under the conditions stated in the theorem the positive operators 
L — QQ^ or L — Q'^Q have spectrum in Ziil?^ such that the point X = is an 
infinitely degenerate point of the discrete spectrum. The eigenf unctions of this 
point satisfy the equation of black or white triangles 

1) QV' = or 2) Q+ij, = 0. 

The eigenf unctions of the ground state A = can be determined by explicit 
formulae, which follow from the construction given above. 

This assertion is a difference analogue of results in , |^ , on the ground 
states of the Pauli operator in a continuous purely magnetic case. 

Problems 1. Prove that the constructed spaces of solutions of the equations 
Qip — and Q^t/j — are complete in £2(Z^) among all solutions, and find 
appropriate orthonormal bases. Explain what kind of cases of solubility in 
£2(Z^) of the equations of black or white triangles exist, apart from the solutions 
found above. Even in the class of operators whose coefficients are exponents of 
linear forms lj(n), this problem has not been solved. 

2. Prove that if the conditions of the theorem are satisfied, then the spectrum 
of the operators L = QQ^ and L = Q^Q has a non-trivial gap A between the 
ground state Aq = and the next level Ai > A > 0. 

3. Find the class of operators Q such that the spectrum of the operators 
L — QQ'^ and L — Q^Q begins at the point A = and is continuous, that is, 
the point A = is the bottom of the continuous spectrum (A = 0). 

The last problem is of special interest. Below we consider algebraic condi- 
tions of the type of commutation of the operators Q, Q"*", which supposedly lead 
to the situation discussed in Problem 3. 

To begin with we consider the difference (purely real) g-analogues of the 
continuous Schrodinger-Landau operator in a homogeneous (constant) magnetic 
field. 

Suppose that the operators Q, have the form (105): 

and depend on the parameters 

c = e^^ , d — e'^^ , Ij = IjiUi + Ij2n2- 

Lemma 16 6 Operators Q, of the form (105) satisfy the relation 

QcdQtd - 1 = <iiQt',d'Qc',d' - 1) (114) 

if the matrix lij and the parameters c, d, c' , d' are related by the following 
equalities: 

c = u^^c', d — u^'^d', q = u^'^ , u ~ e^'^'^ , 

2^11 =2^22 -^12+^21. (115) 
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Remark 12 In the notation 

e'" = M, e'^^ = V 

the operator Q is written in the form 

Qc,d = 1 + cu^'v'^^Ti + d{u^/v)'^'u'^^T2. (116) 

Dependence on the parameters {u,v) in (114) is omitted. Such notation was 
used in p3|] , where these operators were introduced for the first time and their 
eigenfunctions were found. 

The proof of the lemma is by direct verification. 

Comparing the results of Lemma 16 and Theorem 4, we obtain the following 
resuh (||). 

Theorem 5 We consider the operators L = QQ^ and L = Q~^Q under the 
conditions (115). The spectrum of these operators in >C2(Z) for < A < 1 is 
purely discrete, infinitely degenerate, and can only lie at the following points: 

a) Xj = l- j > 0, u < 1; (117) 

b) = 1 - , > 0, M>1. (118) 

In the following cases the spectrum of the operator L occupies all points (117), 
and the spectrum of the operator L occupies all points (117) except Aq: 

a') > v^^ > u^^ > 1; a") > max(i;, -w^^) > 1. (119) 

Analogously, the spectrum of the operator L occupies all points (118), and L 
occupies all points except Aq, if one of the following conditions is satisfied: 

b') > t; > It > 1; b") u > max(i;, i;"^) > 1. (120) 

All eigenfunctions discussed in the theorem are obtained from solutions of 
the equations Qip — or Q~^ip = for proper values of the constants c, d, u, v 
by the use of 'creation operators': ^ ^ ^ 

1) Qcdi^Q = 0, Z^rf^o = ,doQco,doV'o = 0, V = Qt^^do ° • ■ ■ ° A_iV'o, 
ico.do^ = (1 - "^'')V'; 

2) Q+i/'o = 0, Lc.dV'o = Qc,dQXd^o = 0, V = Qco,do o • • • o Qck-ud^-ii>o, 

Cj = u^Cj-i, dj = u^dj-i. 

Remark 13 Of the two cases mentioned in the proof of Theorem 4, in the given 
situation only case 1 is realized, in which the eigenfunctions ip^'^'^ that we found 
vanish on the whole subplane. 



42 



The situation described in Theorem 5 is completely analogous to ||2^, that 
is, to Theorem 1 in §3. The difference is that for n = 2 the spaces of solu- 
tions are infinite-dimensional (the completeness of the solutions found on the 
corresponding levels has not been established so far). 

We note that conditions (119), (120) are obtained directly from the condition 
of Theorem 4, but the case c) here is not realized under conditions (115). 

Problem Prove that the spectrum of the operators L = Qc,dQtd ^'^'^ ^ ~ 
Q^jQcd in the space >C2(^^) for A > 1 is continuous and runs over the whole 
semi-axis A > 1. As we assume, this is a simple Lebesgue spectrum. 

We point out that our operators Q depend on the constants c, d, u, v (if 
the dependence on u, v is not written explicitly, then in the given formula we 
assume that u, v have not changed). The following formulae are true: 

TiQcdTi ^ — Qc'.d', 
T2QXdT.^^^Qc"..d", 

where 

Qc,d = l + ce'i(")Ti+de'=(")r2, 

The same is also true for the operators L = QQ^, L = Q^Q: 

TiLcdTi ^ = Lc'.d', 
T2Lc,dT2 ^ = Lc'i.d"- 

Definition 11 A characteristic vector- section is a function tpcdin) that is ana- 
lytically dependent on all the parameters c, d, u, v and such that 

Lip — X-tJj or Lip — Xip, 

where A does not depend on the parameters. A Block vector- section ipcdi^jp) 
is a function ip such that 

Lrp ^ Xip, TiV'cd = e'^'V'c'.d', T2ipc,d = e^'^'^ipc" ,d" ■ 



c' = cd^^\ d' = de''\ 
c" = cd^^\ d" = de^^^, 



(121) 



By (121) the functions TjTpc,d satisfy these equations, but with shifted param- 
eters (c', d') or (c", d"). This definition makes sense also in the one-dimensional 
case (54), where we had 

= 1 + ce'(")T, l{n)^ln, 



T = Ti, n = ni, TiQcT^^^Q, 



ca 7 



a — eK 



Let us construct the Bloch vector-section for A = when a > 1. We choose 
c — cq > 0. Choosing an arbitrary ipo.coj Put 

^coa'^in,p) = Tr{Mn))e-^^P. (122) 
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Limitations on the construction appear as a ^ 1. There arc no sohitions at all 
in £2 (Z) for a < 1 if we are talking about the equation QcV' = corresponding 
to A = 0. Therefore, a > 1, c 7^ 0. 

As a ^ 1 + the operator tends to an operator with constant coefficients, 
where there are singularities of the Bloch vector-section. 

Apparently, as |c| — > 1, a ^ 1 we can talk about the transition of the Bloch 
vector-section corresponding to A < 1 into an ordinary Bloch solution for the 
equation with constant coefficients (that is, a solution of pure exponent type), 
since for c 7^ the spectrum of the operator for a = 1 has the form (1 — |c|)^ < 
A<(l + |c|)2. 

In addition, for a = 1 the spectrum is bounded in £2, |A| < (1 -|- |c|)^, and 
for a > 1 the operator is of course unbounded. Hence, we necessarily have a 
singularity here. For A > 1 the situation is not clear. 

Problem For the operators L = QQ'^ or L = Q^Q., whose coefficients are 
exponents of linear functions, construct a complete basis of eigenfunctions in 
the form of 'Bloch sections' that depend analytically on all the parameters 

c, d, lij, which has minimally possible singularities in the submanifolds, where 
we have operators with constant coefficients in one of the variables. 

We return now to the case when the point A = is supposedly the lower 
bound of the continuous spectrum for operators of the form 

L = QQ+, (123) 

where the Q have the form (105). 

Hypothesis If the operators Q, commute, QQ~^ = Q^Q, and they have 
a general solution Qij) = 0, Q~^4> = that is bounded on the lattice (or 
is growing sufficiently slowly?), then the point A = is the lower bound of 
the continuous spectrum in £2(Z^) for the operator L = L. We state below a 
natural generalization of this hypothesis. 

Example 10 Let the operators Qc,d) Qtd such that In = Z22 = or w = 1. 
Then we have 

as follows from (114). The operators have the form 
Qc,d = l+cv"=Ti+dw-"iT2, 

We pose the following question. Let two operators be given in the form 

Qi = 1 + a„Ti + 6„T2, Q2 = 1 + c„Tf 1 + dnT^'^ (126) 
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with non-vanishing coefficients. We consider the system of equations of black 
and white triangles simultaneously: 



M = (127) 

In this case is the system (127) completely locally consistent? 

Richer in content is the following formulation. Let the 'initial condition' for 
(127) be given in the form of two arbitrary values of ijj at the ends of any edge 
of the lattice (for example, "(An-Ti and ijjn)- Using equations (127), we can solve 
the system and find the value at any other point n' G 1? ^ moving along the 
paths of the triangles, black and white: knowing the values at the ends of any 
edge, we find the value at the third vertex of any black or white triangle that 
has this edge as its border. Considering a bundle of six triangles with a common 
vertex n, we can pass from one cycle and return again to the edge [n — Ti,n], 
and the linear 'transformation of curvature' will be given by an upper triangular 
(2 X 2)-matrix: 

Anew) _ ■ Aold) „ . (old) 
Yn-Ti — ^nWn-Ti + ^nWn > 

The matrix coefficients are the 'curvatures' An, Bn- They are expressed in terms 
of the coefficients of the system (127): 



A-n 



bnCn+T2dn0.n-T2 



Bn = " ( 1~( ""'"'^^ (3 bn-T2 ) — Cn+Ti + 1 1 " C^n+Ta 1 — dn-Ti 

Cn+T2 \bn \an-T2 \dn 



(128) 



Lemma 17 7 The condition ^„ = 1, B„ = for all n €: I? is necessary 

sufficient for the possibility of a unique solution of the system (127) on the plane 
1? under arbitrary initial conditions imposed at the ends of any fixed edge. For 
this it is also necessary and sufficient to satisfy the following algebraic relation 
for the operators Qi, Q2- 

((Qi-1)(Q2-1)-1) =/„((Q2-l)(Qi-l)-l), fnJ^O. (129) 

We omit the proof of the lemma. 

Corollary 2 // the system (127) is completely consistent, then its solution is 
determined by an arbitrary solution of a one- dimensional difference equation 
of the second order along any path of edges without self-crossing, in particular, 
along the path n\ = const {we obtain an equation for the variable n^) or n2 = 
const {we obtain an equation for the variable n\). If the curvature is trivial, 
then all these equations are equivalent. 
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Example 11 For the commuting operators Qi = Q2 — Qc.d from (125) the 
condition for consistency is satisfied. For the variable ni (712 = const) we get 
the equation 

ct"^ ii^n-T, + i^n+T, ) + (1 + c'^''"^ " ^^^"'"^ )V'n = 0. (130) 



To end this section we give some information on factorizations and Laplace 
transformations for general non-self-adjoint operators L of the second order on 
an equilateral triangular lattice, and also on complex Hermite operators. The 
following assertions hold. 

Theorem 6 To represent the operator 

L = an+bnTi+CnT2+dn-T,Tr^T2 + en-T,Tr^+fn-T,T^^+gn-T,TiT2^ (131) 

in the form 

a) L = Q1Q2 +w„ or h) L = Q2Q1 + w„, 
where Qi and Q2 have the form 

Ql=Xn+ VnTl + ZnT2, Q2 = Pn + Qn-TiT^^ + rn-Ts^TS 

it is necessary and sufficient to satisfy the conditions 

a) bn+T2dnfn+Ti = gnen+T2Cn+Ti, (132) 

b) fndnbn — CnCnQn- 

If both conditions are satisfied, then the Laplace transformations 

a) L = QiQ2+w\ — >L = Q2W-^Qi + l, (133) 

b) L = Q2Qi+v^L = Qtw-^Q2 + l (134) 
can be iterated, moreover infinitely many times [if w,v ^ 0). 

Remark 14 The factorization conditions (132) are invariant with respect to 
the rotation of the lattice by the angle 27r/3. By analogy with Definition 10 
wc can introduce six types of Laplace transformations. Then the 1st, 3rd and 
5th transformations are unitarily adjoint and inverse to the 4th, 6th and 2nd 
transformations respectively. The same applies to the following assertion. 

Theorem 7 To factorize an operator L of the form 

L = an + bnTi+CnT2+dn-TiT2T{~^+bn-TiT^^+Cn-T2T2^+dn-T2TlT2^ (135) 

in the form 

a) L = QQ~^+w„ or b) L = Q~^Q + w„, 
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it is necessary and sufficient to satisfy the following conditions: 

a) 6n-TiC„_T2<in-Ti-T2 € M, 

(136) 

b) d„c„bn e M. 

In particular, if both conditions are satisfied, then the Laplace transformations 
can be iterated infinitely many times. 

The proof of these theorems is straightforward. 

For complex Hermite operators L it is natural to define a class of 'phase' 
gauge transformations 

such that /„ G R. These transformations keep the operator formally Hermitian 
with respect to the previous standard scalar product in LiiJ?') 

{^,^)=Y.^^n- (137) 
The operator L naturally has 'real' and 'phase' projections 

L$ = a„ + -B„Ti + CnT-i + Dn-TiT2T^ ^ + Bn-XiT-^ ^ + C'n-T2l2 ^ + Dn-T2TlT2 ^, 

bn = pnBn, Cn = jnCn, dn = SnDn, whcrc /?„ , 7„ , (5„, a„, w„ are real numbers 
and Bn, C'n, Dn are complex with modulus equal to 1. 

Definition 12 A physical magnetic field is the phase part of the product of 
the coefficients of an operator (of connection) along the border of any black or 
white triangle 

e"'*" = Bn-TiCn-T2Dn-Ti-T2, (black) 
e**^' = DnCnBn- (wkite) 



The following lemma is obvious. 

Lemma 18 8 The operator L is reduced by the phase gauge transformation 
e*^^Le~^^^ to a purely real form if and only if the magnetic field (that is, all 
'magnetic fluxes' ^n^) is trivial: 

$«=0, riGZ^, J = 1,2. 



Corollary 3 A complex Hermite operator L admits an unbounded number of 
Laplace transformations of all six forms if and only if it is reduced by a phase 
gauge transformation to a purely real operator. 
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Thus, the g-analogucs of the Schrodinger-Landau operator on the lattice 
have no connection with the magnetic field. Difference operators in a physical 
magnetic field do not factorize, generally speaking. A proper analogue of an 
operator in a homogeneous magnetic field can (having chosen a specific gauge 
transformation, that is, a 'vector-potential' composed of coefficients 6„,c„,d„ 
such that \bn\ = |c„| = \dn\ = 1) be written in the form 

L = 6 - e**"=Ti - e-**"ir2 - e'*("i+"=)T2Tr^ 

§7. Factorizations and Laplace transformations 
on many- dimensional lattices of regular tetahedra 
in 

We consider a lattice in such that the ends of the basis vectors Ti, . . . , Tjv 
of the lattice together with the point form a regular TV-dimensional simplex 
(a tetrahedron for TV = 3). We call this simplex, and also all others obtained 
from it through translations by a vector of the lattice, black tetrahedra. White 
tetrahedra (A^-dimcnsional simplexcs) are the tetrahedron (0, — Ti, . . . , — Tjv) 
and all others obtained from it through integer-valued translations. 

Two black (white) tetrahedra can have no more than one vertex in common. 
A black and a white tetrahedron can have the longest edge in common. This is 
true for any dimension > 2. 

We consider a real self-adjoint operator of the form 

N 

L = an + '^{bk,nTk + bk,n-TkTk^) + {Ckj,n-Ti TkTf ^ + Ckj,n-TuTjT^^). 

k=l k=jLj 

We assume that the coefficients bk,n > Ckj,n are a real connection, and they are 
all non-zero. 

Theorem 8 The operator L admits a representation in the form 

L = QQ~^ + Wn, 

where 

N 
k=l 

if and only if the following condition for the coefficients bk,n, Ckj,n is satisfied: 

bk,n-Tkbj,n-Tj 



Ckj,n-Tk-Tj 



x„ does not depend on k, j. (139) 



Analogously, the condition for factorization L = Q^Q + Wn is the independence 
of the expression bk,nbj,n/ckj,n on k, j. 
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Remark 15 For = 2 wc do not have any conditions. For iV = 3 wc have 
two conditions for any white tetrahedron: the product of the coefficients of 
connection that sit on pairs of opposite (skew) edges of the tetrahedron is the 
same. In all, the tetrahedron has six edges, which make up three skew pairs, 
hence we have two conditions. If we want to have all forms of factorizations, 
then this condition should also be satisfied for black tetrahedra. 

The proof of the theorem, as before, is by direct verification. 

We now consider factored operators, for which the coefficients of the corre- 
sponding operators Q, have the form of exponents of linear functions, as 
for N = 2: 

Q = l + ce'i(")Ti + rfe'=(")r2 + /e'^^^^Ta, (140) 

N 

For such operators the following results are true. 
Theorem 9 The equations 

1) Qtp = {the equation of black triangles) (141) 

or 

2) Q~^il> = {the equation of white triangles) (142) 

certainly have infinite- dimensional spaces of solutions ip G C2i^^) if the follow- 
ing conditions are satisfied: 
' 1) 

(hi h2 hi\ 

(143) 

/i3, (144) 

^33/ 

(^11 ll2 llZ \ 

ll2 I22—I2I+I12 ^23 — ^21 + ^12 I > 0, ^23 — ^21 + ^12 = ^32 — ^31 +^13; 
^13 ^23 — ^21+^12 ^33 — ^31 + ^13/ 

(145) 

2) the same with Uj —Ufi 
and also the conditions obtained from those stated above by an arbitrary permu- 
tation of the indices 1, 2, 3. 







hi 


h2 h3 




la) 




122 123 






\h3 


I23 ^33 


11 




ll2 


h3\ 


12 


^22 — 


hi+l 


12 ^23 J 


13 




I23 


/33/ 



Corollary ^ // the conditions of the theorem are satisfied, the operators L = 
QQ~^ and L = Q^Q have X = as the point of the discrete spectrum that is 
infinitely degenerate. 
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Problem Prove that the eigenfunctions constructed according to the scheme 
of the proof of Theorem 4 give a complete basis for A = 0. Prove that the 
remaining spectrum is separated from A = by a finite gap A > in the space 

Theorem 10 // the relations 

kj + Iji = h (146) 
are satisfied, then the operators Q, satisfy the relations 

Qc4jQt4J - 1 = 9(QJ,d',/'Qc',d',/' - 1), (147) 

where 

q = e^, c' = e-^c, d' = e'^'d, f' = e-^f. 

This theorem generalizes Lemma 16, formulated for N = 2. 

It is essential to note that the conditions (146) contradict the condition of 
Theorem 9, and we cannot explicitly find solutions of the equations 

Qcdji^ = 0, (148) 

Qlaji' = 0, (149) 

that belong to £2(2^) and are necessary to construct the spectrum of the oper- 
ators L, L. 

Remark 16 As before, we can seek solutions of (141) in the form 

where K2{n) is chosen in such a way that the coefficients of the equation for 

do not contain the variable ni. However, we have not succeeded in reducing the 
resulting equation to one variable. After the substitution 

Xn = If ^ V^n2,n3 

we arrive at the two-dimensional equation of triangles 

A„(pn + BniPn+T2 + Cn^Pn+n = 0. (150) 

Unfortunately, the coefficients of this equation are very complicated, and so far 
we do not know how to find solutions of it explicitly. 

Obviously, one equation (141) can always be solved in the half-space of a 
definite direction if the initial conditions are imposed arbitrarily on vertices in 
any plane that has one of the four forms: 

ni = const, or 712 = const, or 713 = const, or m + n2 + = const. (151) 
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Here the equation of a tetrahedron (as of a triangle for = 2) is treated as the 
evolution equation. It is not reversible, as it is for N = 2, since a solution in the 
inverse direction of 'time' is no longer local; its solubility depends on the initial 
condition belonging to a special functional class. 

Now we consider the condition for consistency of the pair of equations (black 
and white) 

Ql'lp = 0, Ql = I + XnTi + ynT2 + ZnTs, 

1 (152) 



Theorem 11 1 The system (152) is completely consistent, that is, its solution 

is determined by an arbitrary solution of some difference equation of the second 
order in one of the planes (151), if the relation 

(1 - Oi)(l - O2) - 1 = /„((i - 02)(1 - Qi) - 1) (153) 

is satisfied, where fn is a non-zero function [in particular, for /„ = 1 this 
equality is transformed into the commutation condition Q1Q2 = QiQi)- The 
corresponding equation of the second order, for example for the plane ni + 112 + 
ns = const, has the form 

(1 - Qi)(l - Q2)V = V-- 
Proof Since equations (152) have the form 

(i-Qi)V = V, (i-02)V' = V', 

we see that the conditions 



(l-Qi)(l-Q2)V' = V' 

and 

(1 - Q2)(l - Qi)V = V", 

which are equations of the second order in the plane ui + n2 + = const, 
should be equivalent. Hence it follows that there is a non-zero function /„ such 
that 

(1 - Qi)(l - Q2) - 1 = /n((l - Q2)(l - Ql) - 1). 

Thus, we have proved the theorem for the plane rii + ri2 + = const. 

Writing our conditions locally, in the star of any vertex, we see that it is 
invariant with respect to rotations that transform the lattice into itself. All 
four directions of the planes mentioned above are equivalent. We have proved 
Theorem 11. □ 
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Example 12 Let us consider the operators (147), where h = 0. In this case 
the operators Q, commute. The corresponding equation in the hyperplane 
ns = const has the form 

= (A + C^e^""^ + £)2e-2«ni + ^^an, (^^^^^ ^ ) 

(154) 

Now we consider a 'vector factorization' for iV = 3, where the operators 
represent vectors (Q"'""), a = 1,2. Hence, the factorization of the operator is 
given in the form 

2 

a) WW L = ^Q"0+«+«;„, 

a=l 

A (155) 

b) BB L = ^g+"0" + w;„, 

c«=l 

c) BW L = QiQ+i + Q+2g2 + y,^, 

where 



Lemma 19 9 Vector factorization of a self-adjoint operator L is always possible 
{although it is not unique). Moreover, we can regard the potential as a constant 

Wn = Wq. 

The vector factorizations of all three forms generate Laplace transformations 
which, however, cannot be iterated. 

The search for eigenfunctions of the ground state can sometimes be carried 
out in one of them: 

a) Q+V = 0, Q+V = for (WW), 

b) Q^^j = 0, Q'^^j = for (BB), (156) 

c) Q+V = 0, = for (BW). 

If the constant wq for the factorization is chosen correctly, then the lowest level 
is obtained from equations (156), a), b), or c). 

It is difficult to find a criterion for the existence of one solution for the 
systems (156). However, as before, we state here a criterion for complete local 
consistency of these systems. 

Theorem 12 2 The system (156) a) (or b) has solutions that are uniquely 
determined by arbitrary initial conditions imposed on the vertices of any straight 
line with direction vector Ti, T2, Ts or TiT^^ , i ^ j , inside the dihedral angle 
in which lie the white {respectively, black) tetrahedra adjoining this straight line 
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along an edge, if and only if the following condition is saMsfied: for any four 
white {respectively, black) tetahedra that are located in angles of a tetrahedron 
twice the size, the eight equations corresponding to them are linearly dependent. 

Proof For dcfiniteness wc consider the system (156) b). For the values of 
at the vertices of every black tetrahedron we have two equations. This means 
that if we know the values of ip at two vertices of a black tetrahedron, we can 
find the values at the other two. Hence we can easily see that if ijj is given on 
all the vertices of some straight line, then our system dictates the extension of 
ip into the interior of the stated sector. Thus, we only need to clarify when this 
extension is possible for any initial conditions. 

We consider a tetrahedron Tn, similar to the black one, with edge of length TV. 
Let V{N) be the number of black tetrahedra that lie in Tn, V{N) = N{N + 
1){N + 2)/6. For the values of tp at the V{N + 1) vertices that lie inside and on 
the border of Tjv we have 2V{N) equations, of which only V{N + 1) — N — 1 
should be linearly independent in order to satisfy the condition for complete 
consistency. But 2V{N) - V{N + 1) + iV + 1 = V{N - 1) is the number of 
double-size tetrahedra T2 C Tn- □ 

The condition for consistency of the system (156) c) was considered in 
Theorem 11. 

§8. Factorizations of operators and Laplace trans- 
formations on two-dimensional surfaces 

Let us consider a two-dimensional manifold without boundary, triangulated in 

such a way that two-dimensional simplexcs can be painted in two colours (black 
and white) so that two triangles that are adjacent along an edge have different 
colours. In this case an even number of triangles should meet at each vertex. 

The metric on the surface is chosen in such a way that all triangles are 
equivalent to an equilateral triangle in the Euclidean plane. Hence it follows 
immediately that at any vertex the total curvature has the form (27r — Nt^Tr/S), 
where TVtr is the number of triangles meeting at a given vertex. If A^tr = 6, then 
the curvature is equal to zero. If A^tr = 4, then the curvature is positive. If 
A'tr > 6, then the curvature is negative. 

Definition 13 A vertex scalar Schrddinger operator is an operator L that acts 
on a function from a vertex according to the formula 



where summation goes over the vertices P' that are nearest to P but do not 
coincide with P. The condition that the operator L is self-adjoint has the form 




(157) 



p' 



bp'.p = bp- pi 



ap e K. 



(158) 
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Definition 14 1) A triangular {black) Schrddinger operator of type I is an 
operator L that acts on functions tpT of the black triangles T in such a way that 

(LV)t = bT:T'1pT' + dT'lpT, (159) 
T' 

where T' arc all the black triangles that have a common vertex with T. 

2) A triangular {black) Schrddinger operator of type II is an operator L that 
acts on ipT according to formula (159), but now T' runs over all black triangles 
that have with T a common neighbouring (that is, adjacent along an edge) white 
triangle. 

Analogously we define white triangular operators of types I and II. In both 
cases the self-adjointness condition with respect to the standard scalar product 

T 

has the form 

bx'-.T — bx-.T'j ax G K- 



Theorem 13 3 Any real self-adjoint vertex operator L admits a unique factor- 
ization of the form 

L = QQ+ + wp, 

where 

{Q~^iP)t = ^ yT:Pi>P + XTtpT, (161) 
P 

and the sum is taken over all vertices P of the black triangle T . For any black tri- 
angle T and two of its vertices P\, P2 we have the equality {compare with (101)): 

yT:PiyT:P2 = bpr.P2- (162) 

We have also an analogous assertion for white triangles. 

Theorem 14 4 A real triangular black self-adjoint operator L of type I admits 
a factorization of the form _ 

L = Q+Q + VT, (163) 

where the operator has the form (161), if and only if for every vertex P the 
matrix Bp = {bx-.T'), where the black triangles T, T' belong to the star of the 
vertex P, has the form 

Bp = Diag+Ap, rkAp = l. (164) 
For multiplicities N^^ = 4, 6 this condition is always satisfied. 
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The operator Q in (163) is defined up to a transformation that does not 
change any product of the form 

VT-.PVT'-.P (T ^ T'). 

In particular, for Nt^ > 6 the coefficients yr-.p are defined uniquely, and for 
= A up to the transformation (T ^ T') 

Vt-.p ^ IJ-pUt-.p, Vt'-.p i-^ fJ-p^VT'-.p, 
where ^ fip is any non-zero function from vertices with Nii{P) = 4. 



Theorem 15 5 A triangular {black) Schrodinger operator of type II always 
admits a factorization of the form 

L = QQ++UT, (165) 

where 

(<3"^V')ti =^yTi:T'tpT + XTilpTi, 
T 

Ti is a white triangle, and the sum is taken over all black triangles T adjacent 
to Ti. In fact, the coefficients yxi-.T correspond to the edges of the triangulation. 
We have the equalities 

yr-.TiVT'-.Ti = br-.T' (166) 

for any three triangles T, Ti, T' sequentially adjacent to each other {black, white, 
black) from the star of a vertex with iVtr > 6; 

Vt-.TiVt'-.Ti + Ut-.tiUt'-.ti = br-.T' (167) 

for black triangles T, T' and white triangles Ti, T{ from the star of a vertex 
with Ntj. = 4. 

The factorization is unique in the neighbourhood of all vertices P such that 

Nti{P) > 6. In the neighbourhood of vertices P with Ntr{P) = 4 the factoriza- 
tion is not unique: every equation of the form (167) can be replaced by two: 

Vt-.TiVt'-.t^ = ^tVt' ' VT-.TiVT'-.Ti = b^T-T' ' (168) 

where br-.T' = b^}rp, + b^}j„ is an arbitrary decomposition of br-.T' into a sum of 
two non-zero terms. The resulting system is uniquely soluble. 

In this theorem the black and white colours can be interchanged. 

Remark 17 If for all vertices of the triangulation we have A^tr = 4, 6, then the 
classes of operators of types I and II coincide. 
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Corollary 5 1)// for all vertices we have Ntr < 6 {that is, the curvature is 
non-negative) , then factorizations of all types are defined. Correspondingly, all 
the Laplace transformations related to them are also defined. For example, for 
a vertex operator the Laplace transformation has the form 

L = QQ+ + w< — >L = Q+w-'^Q + l. (169) 

2) Factorizations of triangular Schrodinger operators L of type II are always 
defined in terms of triangles of the opposite colour. The degree of non-uniqueness 
of these factorizations and the corresponding Laplace transformations is deter- 
mined by the number of vertices P of positive curvature, N^-r{P) = 4. 

Remark 18 For an equilateral triangular lattice in the plane the spaces of 
functions on vertices, on black and on white triangles, can be identified, and then 
these factorizations coincide (up to equivalence (104)) with the factorizations of 
operators on a regular lattice considered in §6. 

Conclusion. Formulation of the problem on cyclic Laplace chains arises in two 
cases. 

Case 1: for N^-^ < 6 (Corollary 5). 

Case 2: for vertices with A^tr = 4 under the conditions of Corollary 5. 

In case 2 we can limit ourself to Laplace transformations for operators of 
type II, which transform functions on black triangles into functions on white 
triangles and vice versa. This is similar to one-dimensional Darboux transfor- 
mations. 

§9. Simplicial connections. Generalizations 

We consider a simplicial complex K. Everywhere we denote a simplex of di- 
mension I by (t'. 

Definition 15 A simplicial connection of type {q,j, k), < j < k, is an equation 

^ C^^+k.^qlpcrl = 0, (170) 
o-<ico-«+'= 

determined by a vector-function which ascribes to every pair of simplexes 
fT«,CT«+'= C K such that a« C c7«+'= a collection 

c^,+fe,^<, = (c«<,+,^^<,) eM", a = l,...,m, (171) 

where 

™ = ^q+k+l - 

We require that the function c^q+k.^g satisfies the conditions for non-degeneracy 
and localization (see below). 
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Such a function defines an operator which transforms the space of nu- 
merical functions of simplexes a' into the space of m-vector-functions of sim- 
plexes 0"^+'^: 

(Q+V)<..+fc = ^c^,+.,^,V-TO e K™. (172) 

Simplicial connection depends only on the zero-space of this operator: 

g+"V = 0, a=l,...,m. (173) 

In fact, this relation is written separately in each simplex 17*+'^. 

Requirement of non-degeneracy. Equation (197) should be such that for any 
subsimplex 0-'+'^ C K an arbitrarily given collection of values 'ipai in sim- 
plexes a'^, forming a g-dimensional skeleton of any simplex cr'+^ C cr'"'"'^, should 
uniquely and consistently determine values of ip on all remaining gr-dimensional 
subsimplexes in C K. Thus, we can arbitrarily define values of tp in the 
g-dimensional skeleton of simplexes C7*+^ C C7*+'^. 

Lemma 20 Suppose we are given a path 7 composed of {q + k)- dimensional 

simplexes (t'"'"'^, (t|"'''^, . . . , where uf^'^ and crf^f intersect exactly along a face 
of dimension q + j for all s. Then the simplicial connection consistently defines 
a solution of (170) along the path 7, starting from arbitrary initial data defined 
in a q-dimensional skeleton of an arbitrary (q + j)- dimensional face of any of 
the simplexes crf"'"'^ that make up this path. 

We obtain the proof of the lemma in a simple way from the definitions: 

a solution is constructed by transition from the simplex crf"'"'^ to the simplex 

— 10 
(Js+i , s — . . . . 

The analogue of curvature arises naturally for 'closed' paths 7 in which 
(T^"*"*^ = cr^'^. In this case there arises a 'holonomy transformation' of the 
simplicial connection (170) along the path 7: solving (170) along the path ac- 
cording to the lemma, starting from some face ctq"*"-' C o'I^''' , we finally find the 
value of tjj on all g-dimensional faces of the simplex cr^*^, including C cTq^'' ■ 
These values may not coincide with the initial values. There arises the linear 
transformation 



oM , TreM 



where M — C^q+j+i the number of g-dimensional faces in the simplex a^^'' C 
af'^'^ on which the values of the function ipc^q were given arbitrarily. 

Starting from the requirement of non-degeneracy we can reduce equation 
(170) in any simplex 0-'+'^ to such a form that inverse operators are defined for 
any pair of simplexes a\'^\a2^-' C (7^+'^: 

where the space consists of all possible values of on simplexes cr^ C (t|+-' , 
e = 1, 2, and the bases in them are ^- functions of the simplexes cr^. 
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Requirement of localization. All operators L12 are uniquely defined only by the 
minimal simplex containing af^-' and ct^"'"'' : 

This means that if two simplcxcs ap''' and (t|^'^ intersect along the simplex 
g-g+j+s^ where s > 0, then for them all the operators L12 inside the simplex 
g-q+j+s coincide. 

The requirment of localization is automatically satisfied in two cases. 

Case 1. k = j + l. Here always cr«+J+i = cr«+'=. 

Case 2. In a simplicial complex K the simplexes af'^'' and a^'^'^ coincide if they 
intersect in a face of dimension greater than q + j. 

In the examples already considered in this work (above) we have always had 
q = 0, that is, functions ^ were defined on vertices. Localization corresponded 
to cases 1 or 2. 

A trivial example Let k = l,q = Q,j = 0. We have the usual Abelian 
connection, sitting on edges of the complex K. In fact, from (170) 

we arrive in this case at the operator of multiplication 

Ll2 = — V<T? I ' = Ll2{lpa°)- 

Example 13 a) Let g = 0,fc = 2,j = l and let K he & triangulation of a two- 
dimensional surface. Here we have one equation. This situation was studied 

above (sec §6). 

For K = ^ with an equilateral triangular lattice we paid particular at- 
tention to the case of 'zero local curvature' where a solution of the equation 
<5+V = reduces to a difference equation of the second order on the straight 
line (see Example 11). According to a hypothesis of the authors, this situa- 
tion arises in some interesting cases when the point A = is supposedly the 
lower bound of the continuous spectrum of the Schrodinger operator (123). On 
surfaces with non-trivial topology the case of zero local curvature leads to the 
global monodromy defined on the group tti . 

b) Let g = 0, .) = 1, A; = 2 and let K be the 'black part' (or 'white part') of 
the black and white triangulation of a surface. Here we also have one equation 
Q'^ip = 0, but we do not have local curvature, since black triangles are adjacent 
to each other at vertices (there are no common edges). Therefore there are no 
non-trivial paths. 

c) Let g = 0, j = 0, fc = 2, and let K be the same as in b). Here we have 
curvature. 

Example 14 Let g' = 0, A; = 3, and let K be the lattice of regular tetrahedra 
in B?. We have K = K^yj K2, where Ki is the 'black' part, K2 is the 'white' 
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part, and i^i nif2 is a onc-dimcnsional skeleton. Adjacency of two tetrahedra in 
K is possible along no more than an edge, and adjacency of tetrahedra inside Ki 
(or K2) is only possible at a vertex. We have already considered the following 
cases a), b), sec §7: 

a) j = 2, complex Ki (or K2). There is no local curvature here. 

b) j = 2, complex K. There is also no local curvature here, since there are 
no paths 7. Nevertheless, there is a non-trivial condition for consistency (153) 
which ensures an extensive space of solutions of the equation Q^tp = 0. This 
particular analogue of curvature has not been studied in a general form. 

c) j = 1, complex K. Here the concept of local curvature arises, since 
we have many paths 7. The system (170) has a solution in the whole of Z^, 
which is uniquely determined by an arbitrary pair of values of t/j at the vertices 
of any edge if and only if the following condition is satisfied: eight equations, 
corresponding to each set of four tetrahedra, a pair of black and a pair of white, 
from the star of one vertex, such that both black tetrahedra adjoin each of the 
white ones along an edge, are linearly dependent. This assertion means that our 
system in the star of every vertex has a two-dimensional local space of solutions. 
Geometrically, 'curvatures' that obstruct this solution correspond on the border 
of the star — a sphere S*^ divided into triangles and squares just to squares. In 
all we have six conditions at each vertex (the number of squares), but one of 
them is dependent. 

d) Let q = 0, j = 0, and let K he Ki. There arc paths here composed of 
tetrahedra, linked along vertices. The space M*^ for j = is one-dimensional, 
and we arrive only at Abelian connections. This is true for j = in any simplicial 
complexes for any q. 



Lemma 21 1 A simplicial connection of type {q,j,k) that satisfies the con- 
ditions for non- degeneracy and localization defines a multiplicative curvature 
transformation 

Ftj — > M'^, Rry^ry^ = Ftj^Rj^^ Rj—i = R^^ (174) 
for any path 7, 

q+k q-\-k q+k 

where cr^*^ = af^'' and the intersection crf"'"'^ fl cr^^i is a face of dimension 
greater than or equal to q + j for all s = 1, . . . , N . 

The proof of the lemma is easily obtained from the previous one. 
We now consider real discrete self-adjoint operators 

(LV-).! =E^-?-|V'a|, (175) 



acting on functions of q-simplexes of the complex K. We consider only 'operators 
of the second order', where the coefficients b^^.„^ are different from zero only 
for 'nearest neighbours'. 
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Definition 16 The simplcxcs af, (7| are fc+-nearcst (fc_-ncarcst) if this pair is 
contained in some simplex a'^+''+ (respectively, if the intersection fl cr| is a 
simplex ct'"'^-). 

We can represent the operator L as the sum of operators of two types L = 
L+ + L- . We consider these types separately; 

= X] Kl.a-'+''+:al^<T'l' {type+) 

where cr^jCrf are fc+-nearcst and (^1,(^1 are fc_-nearest. Let k+=k- =fc. We 
consider an operator (5+ of the same type as in the definition of a simplicial 
connection, and the operator Qi: 

Here the images are vector-functions, and the "0 are scalars. 

Definition 17 1) A factorization of the first type is a representation of L+ in 
the form 

L+ = QQ++w, (176) 

where (5+ is the operator adjoint to Q and w = w„q is the operator of multipli- 
cation by a numerical function. 

A special factorization of the first type is a representation (176), where w = 
const. 

2) A factorization of the second type is a representation of the operator L_ 
in the form 

L-=Q+Qi+v, (177) 

where v is the operator of multiplication by a function of cr' and the operator 
Qf is adjoint to Qi. 

A special factorization of the second type is the case when v = const. 

3) A Laplace transformation of the operators L± is defined by the formulae 

L+^L+ = Q+w-^Q + l, ^ ^ 

178 

L_=Qiv-^Q+ + 1, 

and their zero eigenvectors are transformed by well-known formulae into the 
eigenvectors of L±: 

^+ I — ,ip+ = Z+^+ = 0, = 0, 
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Obviously, the operator is written analogously to where = v~^/'^Q'^ . 
These formulae assume that w > 0, v > 0, although this does not play a serious 
role here. 

Detailed studies have been devoted above to the investigation of factoriza- 
tions and Laplace transformations in particular cases. Here we consider only 
the simplest example, that is, a standard discrete real Schrodinger operator. 

Example 15 Let q = 0, A: = 1, j = 0. We have an operator defined on functions 
of vertices 

L = L+, {Lii))„o= ^ b^o,„i,„oii)^Q+h„od{al,al)ti)^o. 

(The last term is present only for CTj = (Ji] this is shown by the J- function.) 

We can always factorize such an operator in the form (we always suppose 
that h^o.^i.^o ^ 0) 

L = QQ++w, (179) 
where w is some function. This factorization is not unique. 

Question Is there always a 'special factorization', where w = const? It exists 
for the complex K = Z,^, see §3, but in the general case the answer is not clear 
(we assume that b„o > 0). 

Remark 19 Earlier (see §8) in connection with the Laplace transformation on 
surfaces with black-white triangulation we considered yet another form of the 
factorization: let the complex K and two subcomplexcs Ki, K2 be given, where 
Ki n i4r2 is a (g — l)-dimensional skeleton of the complex K, and the operator 
L acts on functions of cr^ C K\ (white g-simplexes) : 

The notion of the 'nearest fc-neighbourhood' is defined for simplexes from Ki by 
means of faces of dimension q — k, which arc common for them with an arbitrary 
'black' g-simplex from K2. The factorization is sought in the form (179), where 
w is the operator of multiplication by a function, and is an operator of the 
form 

where C K2, cr'' C Ki, cr*~'^ C a"^ n cr"^. Analogously, the operators L act on 
'black' g-simplexes o"^ C K2 and the factorization is sought in the form 

L = Q+Q + v. 

Previous formulae define Laplace transformations and also special Laplace trans- 
formations if the potentials v, w are constant. 
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Example 16 Let be a g-dimcnsional triangulated manifold with black- 
white coloured g-dimensional simplexes, K = KiU K2, and let Ki fl 1^2 be a 
{q — l)-diniensional skeleton. For fc = 1 we have tr'"^ as a common face of 
exactly two simplexes ct'?,(t'^, that is, black and white. By analogy with the 
case q = 2 (see §8), the condition for factorization depends on the iV(c7'~^)- 
multiplicity of {q — 2)-dimensional simplexes cr*~^ C K, which is equal to the 
number of adjacent {q — l)-dimcnsional (t*~^ D (t*~^. This number is even, 
> 4. For a] C i^i, a| C K2 we have: 

a) = c,,,,,-^^,c,|^,|-^^<„ = at' n at', N{a-i-^) > 6; 

Let us consider a): the number of equations is c;qual to the number of 
{q — 2)-dimensional faces of the simplex a^, that is, / = q{q + l)/2. The number 
of unknowns is m{q + 1), where m is the dimension of the vector c^a.,yq-i.-^q. 
Thus we have the condition m{q + 1) > q{q + l)/2 or m > q/2. The scalar 
factorization (m = 1) needs to satisfy additional conditions for q' > 3. 

Example 17 We again consider K, a ((7 + 2)-dimensional manifold with 'black- 
white' coloured (g + 2)-simplexes, K = Ki[jK2, and ifinA'2, a (g-l- l)-skeleton. 
Suppose we are given an operator of some other type, acting on functions of 
g-simplexes from K 

and we seek a 'white' factorization (179), where w is the operator of multipli- 
cation by a function, and 

a" 

We arrive at the relation 

The number of unknowns is equal to m{q + 2)(g + 3)/2 = Im (if c„q+2.^9 is an 
m-vector) and the number of equations is equal to l{l — l)/2. The factorization 
condition thus has the form 

/m>/(Z-l)/2, (181) 

(g + 2){q + 3)/2 = l,m>{l- l)/2. 

For q = l,k = 2 we see that m = 3 or 5 is acceptable. In these cases 
factorization is always possible in the form (179). The possibility of the spe- 
cial factorization w — const is of interest: one would like to clarify this ques- 
tion. If the factorization is special, then the question of zero modes of the 
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operator L = QQ'^ is of interest, that is, the question of the sohibihty in the 
space -0 S £2 of the equation Q'^^' = 0- For m = 3 and m = 5 equations 
of this type determine the connections (if the conditions for non-degeneracy 
and localization arc satisfied). Paths along which 'parallel transport' is re- 
alized and curvature is defined consist of 3-simplexes {q — 1) adjoining each 
other along edges (m = 5) or along faces (m = 3). The requirement of 
localization is automatically satisfied for the case m = 5, K = Ki (white 
tetrahedra). 

If m = 3, then purely white paths do not exist. The curvature is defined by 
the pair of equations: Qfip = (white part), Q2'4' = (black part), that is, 

K = KiU K2. 

Now suppose we are given an arbitrary simplicial connection on a complex 
K of type {q,j,k). 

Definition 18 The local curvature of the vertex d K is the curvature of this 

connection on special paths of the subcomplex Kp^ which is the simplicial star 
of the vertex cr°: we take all simplexes with vertex cr", a' C JsTp, paths of 
the form 

7=(a?+^<7r^...,a«+'==<7?+'=), (182) 

such that = CTV«-l+^ 

A transformation determined only by paths 7 C Kp (182) is called the 
local curvature of the simplicial connection at the point K. 

Hypothesis For manifolds K the curvature is trivial if the local curvatures of 
all the vertices are trivial (we note that in a number of cases above K was not 
a manifold). 
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